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We review the theory of vortices in trapped dilute Bose-Einstein condensates and compare 
theoretical predictions with existing experiments. Mean-field theory based on the time-dependent 
Gross-Pitaevskii equation describes the main features of the vortex states, and its predictions agree 
well with available experimental results. We discuss various properties of a single vortex, including 
its structure, energy, dynamics, normal modes and stability, as well as vortex arrays. When the 
nonuniform condensate contains a vortex, the excitation spectrum includes unstable ( "anomalous" ) 
mode(s) with negative frequency. Trap rotation shifts the normal-mode frequencies and can stabilize 
the vortex. We consider the effect of thermal quasiparticles on vortex normal modes as well as 
possible mechanisms for vortex dissipation. Vortex states in mixtures and spinor condensates are 
also discussed. 

PACS numbers: 03.75.Fi, 03.65.-w, 05.30. Jp, 67.40.Db 
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I. INTRODUCTION 



The recent dramatic achievement of Bose-Einstein condensation in trapped alkali-metal gases at ultra-low tem- 
peratures has stimulated intense experimental and theoretical activity. The atomic Bose-Einstein condensates 
(BECs) differ fundamentally from the helium BEC in several ways. First, BECs in helium are uniform. In contrast, 
the trapping potential that confines an alkali-metal-atom vapor BEC yields a significantly nonuniform density. An- 
other difference is that in bulk superfluid 4 He, measurements of the momentum distribution have shown that the 
low-temperature condensate fraction is ~ 0.1, with the remainder of the particles in finite momentum states [Q,^), 
whereas the low-temperature atomic condensates can be prepared with essentially all atoms in the Bose condensate. 
Finally, the condensates of alkali vapors are pure and dilute (with mean particle density n and n|a| 3 <C 1), so that the 
interactions can be accurately parametrized in terms of a scattering length a (in current experiments, alkali-metal- 
atom BECs are much less dense than air at normal pressure). This situation differs from superfluid 4 He, where the 
relatively high density and strong repulsive interactions greatly complicate the analytical treatments. As a result, 
a relatively simple nonlinear Schrodinger equation (the Gross-Pitaevskii equation) gives a precise description of the 
atomic condensates and their dynamics (at least at low temperatures). One should mention, however, that unlike 
the spinless 4 Hc atoms, alkali atoms have nonzero hyperfine spins, and various forms of spin-gauge effects can be 
important ||. 

Bulk superfluids are distinguished from normal fluids by their ability to support dissipationless flow. Such persistent 
currents are intimately related to the existence of quantized vortices, which are localized phase singularities with an 
integer topological charge. The superfluid vortex is an example of topological defects that are well known in liquid 
helium and in superconductors [^). The occurrence of quantized vortices in superfluids has been the object of 
fundamental theoretical and experimental wor k []lo| - |l4|| . Vortex-like excitations exist in the earth's atmosphere 
in superfluid hadronic matter (neutron stars) |16| , and even in rotating nuclei |]l7f . Examples of other topological 
defects that could exist in dilute gas condensates are "textures" found in Fermi superfluid 3 Hc |isf| , skyrmions []l9 20 
and spin monopoles |2l}| . Vortices in the A and B phases of 3 He are discussed in detail in the review articles p2[23|. 
In superfluid 3 He the Cooper pairs have both orbital and spin angular momentum. These internal quantum numbers 
imply a rich phase diagram of allowed vortex structures, including nonquantized vortices with continuous vorticity 
(see also Refs. |5fl||). 

In the framework of hydrodynamics, the vortices obtained from the Gross-Pitaevskii (GP) equation are analogous 
to vortices in classical fluids |^6). Also the GP equation provides an approximate description of some aspects of 
superfluid behavior of helium, such as the annihilation of vortex rings p7| , the nucleation of vortices p8| , and vortex- 
line reconnection p9| , p0[ . 

The initial studies of trapped Bose condensates concentrated on measuring the energy and condensate fraction, 
along with the lowest-lying collective modes and quantum-mechanical interference effects (see, for example, Ref. pl[). 
Although the possibility of trapped quantized vortices was quickly recognized [ ]32"| , successful experimental verifica- 
tion has taken several years p3| p7| . This review focuses on the behavior of quantized vortices in trapped dilute 
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Bose condensates, emphasizing the qualitative features along with the quantitative comparison between theory and 
experiment. 

The plan of the paper is the following. In Sec. II we discuss the basic formalism of mean-field theory (the time- 
dependent Gross-Pitaevskii equation) that describes dilute Bose-Einstein condensates in the low-temperature limit. 
We summarize properties of vortices in a uniform condensate and also introduce relevant length and energy scales of a 
condensate in a harmonic trap. In Sec. Ill we discuss the structure of stationary vortex states in trapped condensates. 
We analyze the energy of a straight vortex as a function of displacement from the trap center and consider conditions 
of vortex stability when the trap rotates. Also we discuss the recent experimental creation of a single vortex and vortex 
arrays. In Sec. IV we introduce the concept of elementary excitations (the Bogoliubov equations) and analyze the 
lowest (unstable) mode of the vortex for different values of the interaction parameter. We also consider the splitting 
of the condensate normal modes due to presence of a vortex line. 

In Sec. V we investigate the general dynamical behavior of a vortex, based on a time-dependent variational analysis 
and on the method of matched asymptotic expansions. The latter method allows us take into account effects of 
both nonuniform condensate density and vortex curvature. We consider normal modes of a vortex in two- and three- 
dimensional condensates. Also we discuss the energy of a curved vortex line and a nonlinear tilting of a vortex in 
slightly anisotropic condensates. In Sec. VI we analyze the effect of thermal quasiparticles on the vortex normal 
modes and discuss possible mechanisms of vortex dissipation. Also we discuss the influence of vortex generation on 
energy dissipation in superfluids. In Sec. VII we consider vortices in multicomponent condensates and analyze various 
spin-gauge effects. In particular, we focus on the successful method of vortex generation in a two-component system 
that was recently used by the JILA group to create a vortex. In Sec. VIII we draw our conclusions and discuss 
perspectives in the field. 



Bogoliubov's seminal treatment |B8| of a uniform Bose gas at zero temperature emphasized the crucial role of 
(repulsive) interactions both for the structure of the ground state and for the existence of superfluidity. Subsequently, 
Gross [ p9pp| ] and Pitaevskii ElJ] independently considered an inhomogeneous dilute Bose gas, generalizing Bogoliubov's 
approach to include the possibility of nonuniform states, especially quantized vortices. 

An essential feature of a dilute Bose gas at zero temperature is the existence of a macroscopic wave function (an 
"order parameter") 4" that characterizes the Bose condensate. For a uniform system with N particles in a stationary 
box of volume V, the order parameter \& = WNq/V reflects the presence of a macroscopic number Nq of particles in 
the zero-momentum state, with the remaining N' = N — Nq particles distributed among the various excited states 
with k ^ 0. The single-particle states for periodic boundary conditions are plane waves 

y-l/2 e ik-r labeled with the 

wave vector k, and the corresponding creation and annihilation operators and ak obey the usual Bose-Einstein 
commutation relations [a^a^,] = Ok,k'- In the presence of a uniform Bose condensate with k = 0, the ground-state 
expectation value (0^00)0 = No is macroscopic, whereas the ground-state expectation value of the commutator of these 
zero-mode operators ([ao, a ])o necessarily equals 1. Hence the commutator is of order 1/ \/No relative to each separate 
operator, and they can be approximated by classical numbers ao ~ aj ~ \/Nq. This "Bogoliubov" approximation 
identifies these classical fields as the order parameter for the stationary uniform condensate. In contrast, the ground- 
state expectation value for all the other normal modes (aj c ak)o is of order unity, and the associated operators and 
Ok require a full quantum- mechanical treatment. 

The existence of nonuniform states of a dilute Bose gas can be understood by considering a second-quantized 
Hamiltonian 



Here T = -h 2 V 2 /2M is the kinetic energy operator for the particles of mass M, V tI (r) is an external (trap) potential, 
and the interparticle potential has been approximated by a short-range interaction « jj(r-r'), where g is a coupling 
constant with the dimensions of energy x volume. For a dilute cold gas, only binary collisions at low energy are 
relevant, and these collisions are characterized by a single parameter, the s-wave scattering length a, independent 



II. TIME-DEPENDENT GROSS-PITAEVSKII EQUATION 




(1) 




(2) 
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of the details of the two-body potential. An analysis of the scattering by such a potential (see, for example p2| , ^3| ) 
shows that g m 4irah 2 /M. Determinations of the scattering length for the atomic species used in the experiments on 
Bose condensation give: a = 2.75 nm for 23 Na Jj4[ , a = 5.77 nm for 87 Rb J45j] . and a = —1.45 nm for 7 Li [Q. In a 
uniform bulk system, a must be positive to prevent an instability leading to a collapse, but a Bose condensate in an 
external confining trap can remain stable for a < as long as the number of condensed atoms Nq remains below a 
critical value N cr ~ d/\a\, where d is the oscillator length pl]j43| ]. If the interparticle potential is attractive (a < 0), 
the gas tends to increase its density in the trap center to lower the interaction energy. The kinetic energy opposes 
this tendency, and the resulting balance can stabilize inhomogeneous gas. A vortex line located along the trap axis 
reduces the peak central density in the cloud of atoms. Thus a vortex can help stabilize a larger trapped condensate 
with attractive interactions in the sense it can contain a larger number of atoms |4jj ] . 

The time-dependent Heisenberg operator V>(r,i) = exp(iHt/H) ip(r) exp(—iHt/h) obeys the equation of motion 
ihdip(r,t)/dt = [i/>(r,t),-H], which yields a nonlinear operator equation 

ih = ( T + Vti ) ^ t ) + g ft (r, t) ftr, t) ^(r, t). (3) 

The macroscopic occupation of the condensate makes it natural to write the field operator as a sum -0(r, t) = 
\l/(r, t) + (f>(r,t) of a classical field \&(r, t) that characterizes the macroscopic condensate and a quantum field (j>(r,t) 
referring to the remaining noncondensed particles. To leading order, the Bogoliubov approximation omits the quantum 
fluctuations entirely, giving the time-dependent Gross-Pitaevskii (GP) equation f^Utl! 

ih ^&r~ = [ T + ^ r +.9l*(r,t)| 2 ]*(r,i) (4) 

for the condensate wave function ^(r, t). Since ip(r, t) reduces the number of particles by one, its off-diagonal matrix 
element (N — l\tp(r, t)\N) oscillates at a frequency corresponding to the chemical potential [i ~ E (N) — E (N — 1) 
associated with removing one particle from the ground state. Thus the stationary solutions take the form 5'(r,t) = 
^(r) e -^ tt / ?i j where ^(r) obeys the stationary GP equation (frequently identified as a nonlinear Schrodinger equation, 
although the eigenvalue /i is not the energy per particle) 

(T + y tr + 5 |*| 2 )* = ^. (5) 

Apart from very recent work on 85 Rb using a Feshbach resonance to tune a to large positive values p8[ |, essentially 
all studies of trapped atomic gases involve the dilute limit (rl|a| 3 <C 1, where n is the average density of the gas), so 
that depletion of the condensate is small with iV' = N — Nq oc J n\a^N N. Typically n\a\ 3 is always less than 
10~ 3 . Hence most of the particles remain in the condensate, and the difference between the condensate number iVo 
and the total number N can usually be neglected. In this case, the stationary GP equation (||) for the condensate 
wave function follows by minimizing the Hamiltonian functional 

H = JdV [** (T + V tI ) * + i 3 |*| 4 ] , (6) 

subject to a constraint of fixed condensate number Nq — J dV \ ^\ 2 « N (readily included with a Lagrange multiplier 
that is simply the chemical potential (i). 



A. Unbounded Condensate 

The nonlinear Schrodinger equation (H) contains a local self-consistent Hartree potential energy Vff(r) = ,g | \I/ (r) | 2 
arising from the interaction with the other particles at the same point. In an unbounded condensate with V tr = 0, 
the left-hand side of Eq. (||) involves both the kinetic energy T and this repulsive Hartree potential g \^f\ 2 — gn for 
a uniform medium with bulk density n. On dimensional grounds, the balance between these two terms implies a 
"correlation" or "healing" length 

V2Mng y/&^d 

This length characterizes the distance over which the condensate wave function heals back to its bulk value when 
perturbed locally (for example, at a vortex core, where the density vanishes). 
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For a uniform system in a box of volume V, the condensate wave function is ^ = y/NoJV w y/N/V, and 
Eq. (0) shows that the ground-state energy i?o arises solely from the repulsive interparticle energy of the condensate 
Eint ~ 5 gN 2 /V. The bulk chemical potential is then given by 

/ dEn\ ATrah 2 n , , 

/ ' 9n=—^. (8) 



The corresponding pressure follows from the thermodynamic relation 



Finally, the compressibility determines the bulk speed of sound s: 



1 f dp\ gn fi Aitah 2 n 



8 = m{^) = M = M = ~Sf5-' ° r ' e ^™tly, s = ; (10) 

Equations (0) and (|l(]) both indicate that a bulk uniform Bose condensate requires a repulsive interaction (a > 0), 
since otherwise the healing length and the speed of sound become imaginary. 



B. Quantum-Hydrodynamic Description of the Condensate 



It is often instructive to represent the condensate wave function in an equivalent "quantum-hydrodynamic" form 

*(r,t) = |*(r,t)|e iS ( r ' t >, (11) 

with the condensate density 

n(r,i) = |*M)| 2 - (12) 
The corresponding current density j = (7l/2Mi)[\Ii r *V\E r — (V^*)^] automatically assumes a hydrodynamic form 

j(r,i)=n(r,t)v(r,t), (13) 

with an irrotational flow velocity 

v(r,i) = V$(r,i) (14) 

expressed in terms of a velocity potential 

, x HS(r,t) , , 

*(r,t) = — (15) 

Substitute Eq. (|ll]) into the time-dependent GP equation (|4|). The imaginary part yields the familiar continuity 
equation for compressible flow 

Fin 

( 4- + V • (nv) = 0. (16) 
ot 

Correspondingly, the real part constitutes the analog of the Bernoulli equation for this condensate fluid 



Mv 2 + V tI + T^l + gn + M^- = 0. (17) 



1 IU ,,- 



/n 



at 



To interpret this equation, note that the assumption of a zero-temperature condensate implies vanishing entropy; fur- 
thermore, the conventional Bernoulli equation for irrotational compressible isentropic flow can be rewritten as |M9 50 



6 + v 

Mv 2 + U + — - + M— = 0, (18) 



n 
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where U is the external potential energy, e is the energy density and e + p is the enthalpy density. Comparison with 
Eqs. (H) and (0) shows that Eq. ( |l7| ) for the condensate dynamics indeed incorporates the appropriate constitutive 
relations for the enthalpy per particle (e + p)/n — (y/nf^Ty/n + gn. 

As a result, the hydrodynamic form of the time-dependent Gross-Pitaevskii equation in Eqs. ([11]) and ( |l7| ) necessarily 
reproduces all the standard hydrodynamic behavior found for classical irrotational compressible isentropic flow. In 
particular, the dynamics of vortex lines at zero temperature follows from the Kelvin circulation theorem p9| , pp[ , 
namely that each element of the vortex core moves with the local translational velocity induced by all the sources in 
the fluid (self-induced motion for a curved vortex, other vortices, and net applied flow). The only explicitly quantum- 
mechanical feature in Eq. ( |l7| ) is the "quantum kinetic pressure " (y/n)~ 1 Ty/ri] as seen from Eq. (0), this contribution 
determines the healing length £ that will fix the size and structure of the vortex core. 

In classical hydrodynamics, the flow can be considered incompressible when the velocity |i>| is small compared to 
the speed of sound. More generally, classical compressible flow becomes irreversible when the flow becomes supersonic 
because of the emission of sound waves (which are still part of the hydrodynamic formalism). In a dilute Bose gas, 
however, Eqs. ( |l6| ) and (|l7| ) neglect the normal component entirely. As discussed below in Sec. IV. B, the system 
becomes unstable with respect to the emission of quasiparticles once the flow speed exceeds the Landau critical 
velocity (which here is simply the speed of sound). The normal component then plays an essential role and must 
be included in addition to the condensate. In this sense, a dilute Bose gas is intrinsically more complicated than a 
classical compressible fluid. 



C. Vortex Dynamics in Two Dimensions 

Vinen's experiment [ft2| on the dynamics of a long fine wire in rotating superfluid 4 He strikingly confirmed Onsager's 
and Feynman's theoretical prediction of quantized circulation |H],[ll| . These remarkable observations stimulated the 
study of the nonlinear stationary GP equation (|^) in the absence of a confining potential, building on an earlier 
analysis by Ginzburg and Pitaevskii of vortex- like solutions for superfluid 4 He near T\ Gross and Pitaevskii 

independently investigated stationary two-dimensional solutions of the form \&(r) = y / nx(r), where n is the bulk 
density far from the origin. Specifically, they considered axisymmetric solutions 

*(r)=e*/feV (19) 

where (r±,cj>) are two-dimensional cylindrical polar coordinates, and / — ► 1 for r± ^> £. Equations jl^ ) and ( |l5| ) 
immediately give the local circulating flow velocity 

(20) 



Mr i 



which represents circular streamlines with an amplitude that becomes large as r± — ► 0. Comparison of Eqs. © and 
( pp| ) shows that the circulating flow becomes supersonic (v rs s) when r± m £. 

The particular condensate wave function ( |l9| ) describes an infinite straight vortex line with quantized circulation 

« = ^«fl-v=A (21) 

precisely as suggested by Onsager and Feynman po|,pT| . Stokes's theorem then yields h/M = J dS ■ V x v, with the 
corresponding localized vorticity 

Vxv = A^)( ri )i (22) 

Hence the velocity field around a vortex in a dilute Bose condensate is irrotational except for a singularity at the 
origin. 

The kinetic energy per unit length is given by 



2M 2M I 1 1 2M 



dr i / r 2 , 



(23) 
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FIG. 1. Radial wave function f(r±/£) obtained by numerical solution of the stationary GP equation for a straight vortex 
line. 



and the centrifugal barrier in the second term forces the amplitude to vanish linearly within a core of radius « £ (see 
Fig. |l|). This core structure ensures that the particle current density j = nv vanishes and the total kinetic-energy 
density remains finite as rj_ — ► 0. The presence of the vortex produces an additional energy E v per unit length, both 
from the kinetic energy of circulating flow and from the local compression of the fluid. Numerical analysis with the 
GP equation (5lj yields E v w (TrH 2 n/M) In (1.46i?/£ ), where R is an outer cutoff; apart from the additive numerical 
constant, this value is simply the integral of ^Mv 2 n. 

To illustrate that the time-dependent GP equation indeed incorporates the correct classical vortex dynamics, con- 
sider a state of the form 



*(r,i) = V^ q - r x(r-r ) 



ijj,t/H 



(24) 



where \ is the previous stationary solution (|T^) of the GP equation for a quantized vortex, now shifted to the 
instantaneous position ro(t), and u is now a modified chemical potential. The total flow velocity is the sum of a 
uniform velocity vq = Hq/M and the circulating flow around the vortex. Substitute this wave function into the 
time-dependent GP equation (Eh. Since x itself obeys the stationary GP equation (Eh with chemical potential u = gn, 



a straightforward analysis shows that \x 
the condensate. The remaining terms yield 



\Mvl 



gn, where the first term arises from the center of mass motion of 



ih 



dx(r- 



dt 



r ) . u dr 
— - = — in 

dt 



Vx(r - r ) = -ihv ■ Vx(r - r ) 



(25) 



This equation shows that drg (t) /dt — vo , so that the vortex wave function moves rigidly with the applied flow velocity 
vo, correctly reproducing classical irrotational hydrodynamics. 

A similar method applies to the self- induced motion of two well-separated vortices at ri and T2 with |ri — r2 1 3> £; 
in this case, 



*(r, t) = Vn x(r - n) x(r - r 2 ) 



(26) 



represents an approximate solution with fi — ng because there is no net flow velocity at infinity. The density 
n \ f(r — ri)| 2 |/(r — r2)| 2 is essentially constant except near the two vortex cores, and the phase is the sum S(r — r\) + 
S(r — T2) of the two azimuthal angles for the variable r measured from the local vortex cores. Substitution into the 
time-dependent GP equation readily shows that each vortex moves with the velocity induced by the other, for example 

dt M K 2 <Mr= ri 

This method also describes the two-dimensional motion of many well-separated line vortices [|52[|5 
of the many- vortex case in 2D was also studied in 1 54 56| . 



(27) 

The dynamics 
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D. Trapped Condensate 



The usual condition for a uniform dilute gas requires that the interparticle spacing ~ n -1 / 3 be large compared to 
the scattering length (n -1 / 3 >aor na 3 <C 1). The situation is more complicated in the case of a dilute trapped gas, 
because of the three-dimensional harmonic trapping potential Vtr = \M (lo 2 x 2 + w^y 2 + uj 2 z 2 ). The stationary GP 
equation (|s|) provides a convenient approach to study the structure of the condensate in such a harmonic confining 
potential. 

For an ideal noninteracting gas (with g = 0), the states are the familiar harmonic-oscillator wave functions with 
the characteristic spatial scale set by the oscillator lengths dj = y/h/Mivj (j = x, y, and z). In particular, the 
ground-state wave function can be obtained by optimizing the competition between the kinetic energy -E^in = (T) and 
the confining energy Eti = (Vtr), where (• • •) = TV -1 j dV vp* • • • denotes the expectation value for the state with 
the condensate wave function ^S. The situation is more complicated for an interacting system, however, because the 
additional interaction energy E m t = (^sl^l 2 ) provides a new dimensionless parameter. The ratio E- lnt /NhjJo serves 

1/3 

to quantify the effect of the interactions, where ujq — (ui x Lu y uj z ) is the mean oscillator frequency. It is not difficult 
to show that this ratio is of order Na/do for Na/do < 1 where do — y/h/Muo is the mean oscillator length P^ , |3~i| , ^3| , 
and of order (Na/do) 2 / 5 for Na/do 3> 1- Thus the presence of the confining trap significantly alters the physics of the 
problem, for the additional characteristic length do and energy hu>o now imply the existence of two distinct regimes 
of dilute trapped gases: 



1. Near-ideal regime 



In the limit Na/do ^ 1, the condensate states are qualitatively similar to those of an ideal gas in a three-dimensional 
harmonic trap, with ground-state wave function ^(r) oc exp [— | [x 2 /d 2 + y 2 /d 2 + z 2 /d 2 )~\ . The repulsive interactions 
play only a small role, and the condensate dimensions are comparable with the oscillator lengths dj. 



2. Thomas-Fermi regime 

In the opposite limit Na/do 3> 1, which is relevant to current experiments on trapped Bose condensates, the 
repulsive interactions significantly expand the condensate, so that the kinetic energy associated with the density 
variation becomes negligible compared to the trap energy and interaction energy. As a result, the kinetic-energy 
operator T can be omitted in the stationary GP equation ([s]), which yields the Thomas-Fermi (TF) parabolic profile 
for the ground-state density |$2| 

n(r) w |* TF (r)| 2 = [fj,- V tr (r) ] 8 [// — V tr (r) ] = n(0) | 1 - £ |U © ( 1 ~ £ ||] , (28) 

where n(0) = [i/g is the central density and 0(x) denotes the unit positive step function. The resulting ellipsoidal 
three-dimensional density is characterized by two physically different types of parameters: (a) the central density 
n(0) fixed by the chemical potential [note that n(0) plays essentially the same role as the bulk density n does for the 
uniform condensate, where [i = gn], and (b) the three condensate radii 



A? = TrV (29) 
3 Moj 2 



j 

The normalization integral J dVn(r) = N yields the important TF relation |3^] 

N = ^n(O) i? 3 = or, equivalently, § = 15 ^ > 1, (30) 

15 15 drip d 

1 /3 

where i?o = (R x RyRz) is the mean condensate radius. This last equality shows that the repulsive interactions 
expand the mean TF condensate radius Rq proportional to iV 1 / 5 . The TF chemical potential becomes 

H=^Mu^ = ^tiuo% (31) 
"o 
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so that fi ^> huj in this limit. The corresponding ground-state energy E — j^huj^R 2 ,/ dfyN — fM-^ follows 
immediately from the thermodynamic relation fi = dEo/dN . 

The TF limit leads to several important simplifications. For a trapped condensate, it is natural to define the healing 
length (Q) in terms of the central density, with £ = [87m(0) a] -1 / 2 . In the TF limit, this choice implies that 

£ Rq — d 2 , or, equivalently, = <C 1. (32) 

do Rq 

Thus the TF limit provides a clear separation of length scales £ <§; do "C Rq, and the (small) healing length £ 
characterizes the small vortex core. In contrast, the healing length (and vortex-core radius) in the near-ideal limit are 
comparable with do and hence with the size of the condensate. 

The quantum- hydrodynamic equations also simplify in the TF limit, because the quantum kinetic pressure in 
Eq. ( |l7| ) becomes negligible. For the static TF ground-state density given in Eq. (|28|), the small perturbations n' in 
the density and in the velocity potential can be combined to yield the generalized wave equation |57fl 

d 2 rt' d 2 n' 
M^ r = V-[(fi-V tr )Vn'} or, equivalently, = V • [s 2 (r)Vn'] , (33) 

where s 2 (r) = [fi — Vt r (r)] /M defines a spatially varying local sound speed. Stringari has used this equation to analyze 
the low-lying normal modes of the TF condensate, and several experimental studies have verified these predictions in 
considerable detail (see, for example, Ref. PH). 



III. STATIC VORTEX STATES 

In the context of rotating superfluid 4 He, Feynman Jll| noted that solid-body rotation with v sh = x r has 
constant vorticity V x v s b = 2f2. Since each quantized vortex line in rotating superfluid 4 He has an identical localized 
vorticity associated with the singular circulating flow (^2j), he argued that a uniform array of vortices can "mimic" 
solid-body rotation on average, even though the flow is strictly irrotational away from the cores. He then considered 
the circulation T = § c d\ ■ v along a closed contour C enclosing a large number Af v of vortices. The quantization 
of circulation ensures that T — J\f v ■ n, where n — h/M is the quantum of circulation. If the vortex array mimics 
solid-body rotation, however, the circulation should also be T — 20, ■ A v , where A v is the area enclosed by the contour 
C . In this way, the areal vortex density in a rotating superfluid becomes 

Af v 2il 

nv = A- = -- (34) 

Equivalently, the area per vortex is simply l/n v — k/2CI, which decreases with increasing rotation speed. Note that 
Eq. ( |34| ) is directly analogous to the density of vortices (flux lines) n v = i?/$o m a type-II superconductor, where B 
is the magnetic flux density and $o = h/2e is the quantum of magnetic flux in SI units (see, for example, Ref. f58f) 



A. Structure of Single Trapped Vortex 

1. Axisymmetric trap 

Consider an axisymmetric trap with oscillator frequencies to z and u>± and axial asymmetry parameter A = lu z /uj± 
(note that A < 1 yields an elongated cigar-shape condensate, and A > 1 yields a flattened disk-shape condensate). 
The conservation of angular momentum L z allows a simple classification of the states of the condensate. For example, 
the macroscopic wave function for a singly quantized vortex located along the z-axis takes the form 

¥(r) =e^|*(r_L,z)|- (35) 

The circulating velocity is identical to Eq. (|2^) , and the centrifugal energy [compare Eq. (]23|)] gives rise to an additional 
term ^Mv 2 = h 2 /2Mr 2 L in the GP equation (|^). In principle, a g-fold vortex with ^ cx e lq ^ also satisfies the GP 
equation, but the corresponding energy increases like q 2 [compare the discussion below Eq. (p3|)]; consequently, a 
multiply quantized vortex is expected to be unstable with respect to the formation of q singly quantized vortices. 
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FIG. 2. Contour plot in the xz plane for a condensate with 10 4 87 Rb atoms containing a vortex along the z axis. The 
trap is spherical and distances are in units of the oscillator length d = 0.791 fim. The interaction parameter is Na/d = 72.3. 
Luminosity is proportional to density, the white area being the most dense. 
(Taken from Ref. 



For a noninteracting gas in an axisymmetric trap, the condensate wave function for a singly quantized vortex on the 
symmetry axis involves the first excited radial harmonic-oscillator state with the noninteracting condensate vortex 
wave function 



^(r) cx e 1 ^ rj_ exp 



rf 2 



(36) 



of the anticipated form ( |35| ) . The inclusion of interactions for a singly quantized vortex in small to medium axisym- 
metric condensates with Na/do < 1 requires numerical analysis ^,^9|. Some phases of rotating BEC in a spherically 
symmetric harmonic well in the near- ideal-gas limit (£ > do) were considered by Wilkin and Gunn |p0[| . By exact 
calculation of wave functions and energies for small number of particles, they show that the ground state in a rotat- 
ing trap is reminiscent of those found in the fractional quantum Hall effect. These states include "condensates" of 
composite bosons of the atoms attached to an integer number of quanta of angular momenta, as well as the Laughlin 
and Pfaffian |^l|] states. 

In general, the density for a central vortex vanishes along the symmetry axis, and the core radius increases away 
from the center of the trap, yielding a toroidal condensate density (see Fig. ^J). This behavior is particularly evident 
for a vortex in the TF limit Na/do ^ 1, when 



n(r u z)^n(0) 1- 



e 



z 

m 



e i- 



e 



(37) 



Here, the density differs from Eq. (|2^) for an axisymmetric vortex-free TF condensate only because of the dimensionless 
centrifugal barrier £ 2 / r l- This term forces the density to vanish within a core whose characteristic radius is £ in the 
equatorial region \z\ <C R z and then flares out with increasing \z\. The TF separation of length scales ensures that the 
vortex affects the density only the immediate vicinity of the core [^7 62 6^] ; this behavior can usually be approximated 
with a short-distance cutoff. For such a quantized TF vortex, the chemical potential /ii differs from /io for a vortex-free 
TF condensate by small fractional corrections of order (do/Ro) 4 In (Ro/do). 



2. Nonaxisymmetric trap 

If a singly quantized vortex is oriented along the z axis of a nonaxisymmetric trap (R x ^ R y ) the condensate wave 
function is no longer an eigenfunction of the angular momentum operator L z . In the TF limit near the trap center 
the phase S of the condensate wave function has the form p4j: 

and the condensate velocity is 



S^-^-^iln^lsin^), (38) 
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v *4 { 7i - \{m - ln (fe) H 2 ^ + sin ( 2 ^] } < < 39 ) 

where = 2R x Ry/ (R x + Ry). Near the vortex core the condensate wave function and the condensate velocity 
possess cylindrical symmetry, while far from the vortex core the condensate velocity adjusts to the anisotropy of the 
trap and becomes asymmetric. 



B. Thermodynamic Critical Angular Velocity for Vortex Stability 

If the condensate is in rotational equilibrium at an angular velocity f2 around the z axis, the integrand of the 
GP Hamiltonian (Q) acquires an additional term — ^f*QL z ^ p5| , where L z — xp y — yp x — —ih(xd y — yd x ) is the z 
component of the angular-momentum operator. Thus the Hamiltonian H' in the rotating frame becomes 

H' = H — QL Z = J dV [** (T + Vtr - QL Z ) * + i«?|^| 4 ] , (40) 

where the variables in the integrand are now those in the rotating frame. Similarly, the GP equations (^]) and (j^) 
acquire an additional term — riL z ^>. 



1. Axisymmetric trap 

The situation is especially simple for an axisymmetric trap, where the states can be labeled by the eigenvalues of L z . 
For example, the energy of a vortex-free condensate E' (Q) in the rotating frame is numerically equal to the energy 
Eq in the laboratory frame because the corresponding angular momentum vanishes. A singly quantized vortex along 
the trap axis has the total angular momentum Nh, so that the corresponding energy of the system in the rotating 
frame is E[(fl) = E\ — NM1. The difference between these two energies is the increased energy 

AE'(Q) = E[(Q) ~ E' (Q) = E x -Eq- Nhfl (41) 

associated with the formation of the vortex at an angular velocity SI. In the laboratory frame (Jl = 0), it is clear 
that Ei > Eq because of the added kinetic energy of the circulating flow. If the condensate is in equilibrium in the 
rotating frame, however, E[(£V) decreases linearly with increasing SI, and the relative energy of the vortex vanishes at 
a "thermodynamic" critical angular velocity Sl c determined by A£"(S1 C ) = 0. Equation ( pT| ) immediately yields 

= Ei-Eo 

expressed solely in terms of energy of a condensate with and without the vortex evaluated in the laboratory frame. 

For a noninteracting trapped gas, the difference E\ — Eq = Nhtu± follows immediately from the excitation energy 
for the singly quantized vortex in Eq. ( |36| ) relative to the stationary ground state. In this noninteracting case, Eq. ( p^ ) 
gives f2 c = luj_, so that the noninteracting thermodynamic critical angular velocity is just the radial trap frequency. 
Indeed, the same critical angular velocity value also applies to a g-fold vortex in a noninteracting condensate, be- 
cause of the special form of the noninteracting excitation energy E q — Eq = Nqhuj± and the corresponding angular 
momentum Nqh. Thus the noninteracting condensate becomes massively degenerate as Q — ► luj_ |36|,^]. Physically, 
this degeneracy reflects the cancellation between the centrifugal potential — ^MU 2 r^_ and the radial trap potential 
\MLo\r\ as Q -> u±. 

Numerical analysis pTj for small and medium values of Na/do shows that Q c /ujj_ decreases with increasing TV, 



and a perturbation analysis |67|68|] confirms this behavior for a weakly interacting system, with the analytical result 
fl c /u}± w 1 — l/(2v / 27r) (Na/d z ) for small values of the interaction parameter Na/d z . Figure || shows the behavior of 
Q C (N) in a spherical trap, based on numerical analysis of the GP equation with parameters relevant for 87 Rb pTf . 

In the strongly interacting (TF) limit, the chemical potential Hi(N) for a condensate containing a singly quantized 
vortex can be evaluated with Eq. (|37|), and the thermodynamic identity Hi = dEi/dN then yields Ei(N). Use of the 



corresponding expressions for the vortex-free condensate gives the approximate expression ]6^ , |69| , [70[ 

5 { 67-/? \ 

fl c w - In I ^ I for a TF condensate. (43) 

2 MR i V s / 
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FIG. 3. Thermodynamic critical angular velocity Q c 
d = 0.791 ^m and N atoms of 87 Rb. 
(Taken from Ref. 



for the formation of a singly quantized vortex in a spherical trap with 



This expression exceeds the usual estimate |1J] Sl c w (h/MR\) ln(1.46-Rj_/^) for uniform superfluid in a rotating 
cylinder of radius R± because the nonuniform density in the trapped gas reduces the total angular momentum 
relative to that for a uniform fluid. Equation ( fl3| ) has the equivalent form 



LUJ_ 



5 d\ , /0.67i?_L\ 



(44) 



This ratio is small in the TF limit, because d\jR\ ~ £/R± <C 1. For an axisymmetric condensate with axial 
asymmetry A = uj z /uj±, the TF relation d\jR\ = (d±/15Na\) 2 ^ 5 shows how this ratio scales with N and A. 

In contrast to the case for repulsive interactions, the thermodynamic critical angular velocity Q c for the vortex 
state with attractive interactions increases as the number of atoms grows |7l],^7j. Since f2 c = u>± for a noninteracting 
condensate, fi c for a vortex in a condensate with attractive interactions necessarily exceeds luj_. The stability or 
metastability of such a vortex is unclear because Q = lu± is also the limit of mechanical stability for a noninteracting 
condensate. 

Approximately the same functional relationship holds between the thermodynamic critical frequency £l c and the 
number of atoms in the condensate Nq Jt2[ for nonzero temperatures. A new feature, however, is that the number of 
atoms in the condensate becomes temperature-dependent: 



iVo 
N 



(45) 



where T c is the critical temperature of Bose condensation. If the trap rotates at an angular velocity f2, the distribution 
function of the thermal atoms changes due to the centrifugal force. As a result the critical temperature decreases 
according to [f72| 



TP 



= 1- 



2 



1/3 



(46) 



where T° is the critical temperature in the absence of rotation. Equations (44)-(^6[) allows one to calculate the critical 
temperature T v (£l), below which the vortex corresponds to a stable configuration in a trap rotating with frequency tt. 
In Fig. [| we show the critical curves T C (Q) and T v (fl). For temperatures below T C (Q) the gas exhibits Bose-Einstein 
condensation. Only for temperatures below T V (Q,) does the vortex state become thermodynamically stable. From 
Fig. ^ one can see that the critical temperature for the creation of stable vortices exhibits a maximum as a function 
of O. 
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FIG. 4. Phase diagram for vortices in a harmonically trapped Bose gas, TV = 10 4 , a/d± — 7.36 x 10~ 3 and A = 1. 
(Taken from Ref. @). 



2. Nonaxisymmetric trap 

A rotating nonaxisymmetric trap introduces significant new physics, because the moving walls induce an irrotational 



flow velocity even in the absence of a vortex 1 49 7$-|76f| . In the simplest case of a classical uniform fluid in a rotating 



elliptical cylinder, the instantaneous induced velocity potential in the laboratory frame is |4S||73j 

A 2 - B 2 

® A = Qxy JpT&> (47) 

where A and B are the semi-axes of the elliptical cylinder. The induced angular momentum and kinetic energy are 
reduced from the usual solid-body values by the factor Io/I s b = [(A 2 — B 2 )/(A 2 + B 2 )\ . In the extreme case B <C A, 
the moment of inertia can approach the solid-body value, even though the flow is everywhere irrotational. 

The thermodynamic critical angular velocity Q c for vortex creation in the same uniform classical fluid depends on 
the asymmetry ratio B/A jr3], and experiments on superfluid 4 He confirm the theoretical predictions in considerable 
detail |77|| . In the limit B<i, a detailed calculation shows that £l c (h/2MB 2 ) ln(-B/£); the appearance of B here 
is readily understood from Feynman's picture of a vortex occupying an area » h/2MQ, [compare Eq. ([mJ)] and hence 
having to fit the area irB 2 fixed by the smaller lateral dimension B. 

The preceding analysis for an axisymmetric dilute trapped Bose gas can be generalized to treat the TF limit in 
a totally anisotropic disk-shape harmonic trap with oj 2 + lj 2 -C lo 2 , starting from Eq. ( fiof ) for the Hamiltonian in 
the rotating frame |7S| ]. The presence of a vortex leaves the TF condensate density essentially unchanged, and this 
Hamiltonian can serve as an energy functional to determine the phase S and hence the superfluid motion of the 
condensate. Since R x , R y 3> R z , the curvature of the vortex is negligible. Hence we consider a singly quantized 
straight vortex displaced laterally from the center of the rotating trap to a transverse position ro = (xo,yo) that 
serves as a new origin of coordinates. The condensate wave function then has the form 

* = |*|e^ +iSo , (48) 

where <f> in the first term is the polar angle around the vortex axis and So is a periodic function of <f>. Varying the 
Hamiltonian gives an Euler-Lagrange equation for .So, and it can be well approximated by the solution for a vortex- free 
condensate, which is M/H times the classical expression (^) with A and B replaced by the TF radii R x and R y given 
in Eq. (p9|), and with x and y shifted to the new origin. 

As in Eq. ( ^ll) for an axisymmetric trap, AE'(xo, yoi ^) gives the increased energy in the rotating frame associated 
with the presence of the straight vortex. A detailed calculation with logarithmic accuracy yields p8[ 
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FIG. 5. Energy [in units of Ai7'(0,0)] associated with a singly quantized straight vortex in a rotating asymmetric trap 
in the TF limit as a function of a fractional vortex displacement Co from the symmetry axis. Different curves represent different 
fixed values of the external angular velocity Q: (a) Q = (unstable); (b) = O m [given in Eq. (^)] (onset of metastability 
at the origin); (c) Q = f2 c [given in Eq. (p2|)] (onset of stability at origin); (d) O = §fi c , where the thin barrier hinders vortex 
tunneling from the surface. 
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(49) 



where Co = x\jR 2 x + y 2 /Ry < 1 is a dimensionless displacement of the vortex from the trap center. Here, the mean 
transverse condensate radius R± is given by the arithmetic mean of the inverse squared radii 
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4/x 



(50) 



Figure || shows the behavior of AE' (Co , fi) as a function of Co for various fixed values of fi. Curve (a) for = shows 
that the corresponding energy A_E'(£o, O = 0) decreases monotonically with increasing Co, with negative curvature at 
Co = 0. In the absence of dissipation, energy is conserved and the vortex follows an elliptical trajectory at fixed Co 
around the center of the trap along a line Vt r = const. At low but finite temperature, however, the vortex experiences 
weak dissipation; thus it slowly reduces its energy by moving outward along curve (a), executing a spiral trajectory 
in the xy plane. 

With increasing fixed rotation speed O, the function AE'((q,H.) flattens. Curve (b) shows the special case of zero 
curvature at Co = 0. It corresponds to the rotation speed 



3 h , fR± 

7T In — 

2 MR 2 , \ C 



for a disk-shape condensate, 



(51) 



at which angular velocity a central vortex first becomes metastable in a large disk-shape condensate. For f2 < fl m , the 
negative local curvature at Co = means that weak dissipation impels the vortex away from the center. For > Q m , 
however, the positive local curvature means that weak dissipation now impels the vortex back toward the center of 
the trap. In this regime, the central position is locally stable; it is not globally stable, however, because AE'(0, f2) is 
positive for CI w Q m . 

Curve (c) shows that AE'(0, fl c ) vanishes at the thermodynamic critical angular velocity 



2MR\ 
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Rx 



- 3 ^ r , 



for a disk-shape condensate. 



(52) 
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As expected, this expression ( |52| ) reduces to Eq. ( f43| ) in the limit of an axisymmetric disk-shape condensate. For 
£1 > Q c , the central vortex is both locally and globally stable relative to the vortex- free state, and the energy barrier 
near the outer surface of the condensate becomes progressively narrower. Curve (d) illustrates this behavior for 
ft = |r2 c . Eventually, the barrier thickness becomes comparable with the thickness of the boundary layer within 
which the TF approximation fail s ]79[|, and it has been suggested that a vortex might then nucleate spontaneously 
through a surface instability [8C|j76| , |8l[ |. For a two-dimensional condensate, a phase diagram for different critical 
velocities of trap rotation vs. the system parameter an z (n z is the area density) is given in |81|. 



C. Experimental Creation of a single vortex 

The first experimental detection of a vortex involved a nearly spherical 87 Rb TF condensate containing two different 
internal (hyperfine) components |33| that tend to separate into immiscible phases. The JILA group in Boulder 
created the vortex through a somewhat intricate coherent process that controlled the interconversion between the two 



components (discussed below in Sec. |VIlj). In essence, the coupled two-component system acts like an SU{2) spin-i 
system whose topology differs from the usual U{1) complex one-component order parameter ^ familiar from superfluid 
4 He (and conventional BCS superconductivity). Apart from the magnitude |^| that is fixed by the temperature in a 
uniform system, a one-component order parameter has only the phase that varies between and 2tt. This topology is 
that of a circle and yields quantized vorticity to ensure that the order parameter is single- valued flQ|Jll|| . In contrast, 
a two-component system has two degrees of freedom in addition to the overall magnitude; its topology is that of a 
sphere and does not require quantized vorticity. The qualitative difference between the two cases can be understood 
as follows: the single degree of freedom of the one-component order parameter is like a rubber band wrapped around 
a cylinder, while the corresponding two degrees of freedom for the two-component order parameter is like a rubber 
band around the equator of a sphere. The former has a given winding number that can be removed only be cutting it 
(ensuring the quantization of circulation), whereas the latter can be removed simply by pulling it to one of the poles 
(so that there is no quantization). 

The JILA group was able to spin up the condensate by coupling the two components. They then turned off 
the coupling, leaving the system with a residual trapped quantized vortex consisting of one circulating component 
surrounding a nonrotating core of the other component, whose size is determined by the relative fraction of the two 
components. By selective tuning, they can image either component nondestructively |}7|]; Fig. || shows the precession 
of the filled vortex core around the trap center. In addition, an interference procedure allowed them to map the 
variation of the cosine of the phase around the vortex, clearly showing the expected sinusoidal variation (Fig. ^) . 

The JILA group has also been able to remove the component filling the core, in which case they obtain a single- 
component vortex j3?j] . This one-component vortex has a small core size and can only be imaged by expanding both 
the condensate and the core, which becomes visible through its reduced density |$2i[7Q|. They first make an image of 
the two-component vortex, next remove the component filling the core, and then make an image of the one-component 
vortex after a variable time delay. In this way, they can measure the precession rate of the one-component empty-core 
vortex and compare it with theoretical predictions |pq| . The data show no tendency for the core to spiral outward, 
suggesting that the thermal damping is negligible on the time scale of ~ 1 s. 

Separately, the ENS group in Paris observed the formation of one and more vortices in a single-component 87 Rb elon- 
gated cigar-shape TF condensate with a weak nonaxisymmetric deformation that rotates about its long axis |m]-^6| . 
In essence, a static cylindrically symmetric magnetic trap is augmented by a nonaxisymmetric attractive dipole po- 
tential created by a stirring laser beam. The combined potential produces a cigar-shape harmonic trap with a slightly 
anisotropic transverse profile. The transverse anisotropy rotates slowly at a rate Q < 200 Hz. In the first experi- 
ments p4j| , the trap was rotated in the normal state and then cooled, with the clear signal of the vortex shown in Fig. || 
(the trap was turned off, allowing the atomic cloud to expand so that the vortex core becomes visible). This order 
was reversed (cool first, then rotate) in a later series of runs J3(|. In both cases, the observed critical angular velocity 
~ 0.7u± for creating the first (central) vortex was roughly 70% higher than the predicted thermodynamic value fi c 
in Eq. (M3). These observations agree qualitatively with the suggestion that a surface instability might nucleate a 
vortex p0[ , |7^ , ^l| . Alternative explanations of this discrepancy involve the bending modes of the vortex (discussed 
below in Sec. IV.D.4 and V.D.2). 
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FIG. 6. (a) Successive images of a condensate with a vortex. The recorded profile of each trapped condensate is fit with a 
smooth TF distribution (b). The vortex core is the dark region within the bright condensate image, (c) The azimuthal angle 
of the core is determined for each image, and plotted vs. time held in the trap. A linear fit to the data gives a precession 
frequency 1.3(1) Hz. 
(Taken from Ref. M). 
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FIG. 7. Cosine of the phase around the vortex, showing the sinusoidal variation expected for the azimuthal angle. 
(Taken from Ref. @). 




FIG. 8. Optical thickness of the expanded clouds in the transverse direction showing the difference between the states (a) 
without and (b) with a vortex. 
(Taken from Ref. @). 
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FIG. 9. Dimensionless angular momentum I per particle vs. dimensionless angular velocity Q/u)±. In the figure 
7 = (2/n) 1 ^ 2 aN/d z . Black lines show stable states and gray lines show metastable states. There are no stable or metastablc 
states in the forbidden ranges I = 0-1 and I = 1-1.70. The rotational symmetry of each branch is indicated. The total an- 
gular momentum diverges as Q approaches the maximum angular velocity lj±. Three-dimensional plots of constant density 
show states with two-fold and six-fold symmetry. Reprinted by permission from Nature 397, 327, (1999), ©1999 Macmillan 
Magazines Ltd. 



D. Vortex Arrays 



Under appropriate stabilization conditions, such as steady applied rotation, vortices can form a regular array. 
In a rotating uniform superfluid, the quantized vortex lines parallel to the axis of rotation form a lattice. This 
lattice rotates as a whole around the axis of rotation, thus simulating rigid rotation flS4j| . At nonzero temperature, 
dissipative mutual friction from the normal component ensures that the array rotates with the same angular velocity 
as the container. Early experiments on rotating superfluid 4 He |85|]l3| , |86|| provided memorable "photographs" of 
vortex lines and arrays with relatively small numbers of vortices, in qualitative agreement with analytical |87 §3] and 



numerical |89 9fJ predictions. A triangular array is favored for vortices near the rotation axis of rapidly rotating 
vessels of superfluid helium [^7j . Vortex lattices also occur in the neutron superfluid in rotating neutron stars |l|] . 

Even before the recent observation of vortex arrays in an elongated rotating trapped condensate |34|,|3^] , several theo- 
retical groups had analyzed many of the expected properties. In a weakly interacting (near-ideal) axisymmetric conden- 
sate, the thermodynamic critical angular velocity Q c for the appearance of the first vortex is already close to the radial 
trap frequency lu±, so that the creation of additional vortices involves many states 4> m (r±) cx e lm ^ r™ exp(— \r\ld\) 
with low energy mh(uj± — Q) per particle in the rotating frame. Butts and Rokhsar |6?]] used a linear combination of 
these nearly degenerate states as a variational condensate wave function, minimizing the total energy in the laboratory 
frame E\ a h subject to the condition of fixed number N of particles and fixed angular momentum I per particle. As ex- 
pected from the theoretical and experimental results for liquid helium, the system undergoes a sequence of transitions 
between states that break rotational symmetry. Several of these have p-fold symmetry where p is a small integer. 
Each vortex represents a node in the condensate wave function, and their positions can vary with the specified angular 
momentum. Indeed, as I increases from to 1, the first vortex moves continuously from the edge of the condensate 
to the center. For larger number of vortices, the centrifugal forces tend to flatten and expand the condensate in the 
radial direction. In this approach of keeping I fixed, the angular velocity follows from the relation Ml = dE^/dl. 
Figure || shows the angular momentum versus the angular velocity for the first several states. Reference |Fl| has 
carried out more detailed studies of the states for relatively small values of the angular momentum per particle I < 2. 

These analyses work at fixed angular momentum Nl, in which case the angular velocity Q must be determined 
from the resulting £'iab(0- I n contrast, the ENS experiments fix f2 (as do experiments on superfluid helium) and then 
measure L z from the splitting of the quadrupole modes ||(| (see Sec. IV. D. 3). The JILA group |92| also uses this 
technique to detect the presence of a vortex in a nonrotating condensate. The transition from fixed L z to fixed f2 can 
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FIG. 10. Arrays of vortices in a Bose-Einstein condensate stirred with a laser beam. 
(Taken from Ref. ||). 



be considered a Legendre transformation to the Hamiltonian (^) m the rotating frame. Even though it is easier to 
work at fixed Q (because there is no constraint of fixed L z /N = I), no such analysis has yet been carried out in the 
weak-coupling limit. 

In the strong-coupling (TF) limit, Castin and Dum ]70| have performed extensive numerical studies of equilibrium 
vortex arrays in two and three dimensions, based on the Hamiltonian in the rotating frame (thus working at fixed f2). 
They also propose an intuitive variational calculation based on a factorization approximation that is very similar to 
Eq. (prj|), apart from a different analytic form of the radial function | B2[ , p3[ . 

The nucleation of vortices and the resulting structures of vortex arrays in zero temperature BECs are also inves- 
tigated numerically by Feder, Clark and Schneider |7(|. In their simulations, vortices are generated by rotating a 
three-dimensional, nonaxisymmetric harmonic trap. Vortices first appear at a rotation frequency significantly larger 
than the critical frequency for vortex stabilization. At higher frequencies, the trap geometry strongly influences the 
structure of the vortex arrays, but the lattices approach triangular arrays at large vortex densities. 

The ENS experiments ]34|,[35| have produced remarkable images of vortex arrays. Figure |l0| shows three different 
arrays with up to 11 vortices (obtained after an expansion of 27 ms). The initial condensate is very elongated (along 
with the vortices), so that the radial expansion predominates once the trap is turned off. As a result, the expanded 
condensate acquires a pancake shape similar to that in Fig. |^. 



IV. BOGOLIUBOV EQUATIONS: STABILITY OF SMALL-AMPLITUDE PERTURBATIONS 

This section considers only the behavior of a dilute one-component Bose gas, for which the analysis of the eigenfre- 
quencies is particularly direct. In the more general case of two interpenetrating species, even a uniform system can 
have imaginary frequencies for sufficiently strong interspecies repulsion [ p3| , p4[ ; this dynamical instability signals the 
onset of phase separation. 



A. General features for nonuniform condensate 



The special character of an elementary excitation in a dilute Bose gas largely arises from the role of the Bose 
condensate that acts as a particle reservoir. This situation is especially familiar in the uniform system, where an 
elementary excitation with wave vector k can arise from the interacting ground state 'J'o either through the creation 
operator a t or through the annihilation operator a_ k (in the thermodynamic limit, these two states ajjj'I'o arid a-k^o 
differ only by a normalization factor). The true excited eigenstates are linear combinations of the two states, and the 
corresponding operator for the Bogoliubov quasiparticle is a weighted linear combination ]3§| , |42] , f43| ] 

a k = Mfe a k + v k a_ k , (53) 

where Uk and Vk are the (real) Bogoliubov coherence factors. This linear transformation ( [53"| ) is canonical if the 
quasiparticle operators also obey Bose-Einstein commutation relations, which readily yields the condition 

u 2 k -vl = l, for all k 0. (54) 

More generally, the second-quantized Bose field operator -ip in Eq. (g) can be written as -0(r) « ^(r) + (j>(r), where 
is a small deviation operator from the macroscopic condensate wave function ty. These deviation operators obey 
the approximate Bose-Einstein commutation relations 
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Since ?A( r ) does not conserve particle number, it is convenient to use a grand canonical ensemble, with the new 
Hamiltonian operator K = H — fj,N instead of the Hamiltonian ([[]) . To leading (second) order in the small deviations, 
the perturbation in K contains not only the usual "diagonal" terms involving but also "off-diagonal" terms 
proportional to 4><j) and (fyefy. Consequently, the resulting Heisenberg operators <j> and c/y obey coupled linear equations 
of motion (it is here that the role of the condensate is evident, for this coupling vanishes if '5 vanishes). Pitaevskii 
developed this approach for the particular case of a vortex line in unbounded condensate, and the formalism was 
subsequently extended to include a general nonuniform condensate 95 9^]. 

In direct analogy to the Bogoliubov transformation for the uniform system, assume the existence of a linear trans- 
formation to quasiparticle operators oij and or- for a set of normal modes labeled by j 



(r,t)=Yl [u 3 {v)a 3 {t) - v*{v)a](t) 



(56a) 



(56b) 



where the primed sum means to omit the condensate mode. Here, the quasiparticle operators a., and a\. obey Bose- 



Einstein commutation relations 



= 5jk and have simple harmonic time dependences otj(t) = otj exp (—iEjt/h) 

shows that the corresponding 



and <Xj(t) — ctj exp (iEjt/H). Comparison with the equations of motion for (j> and 
spatial amplitudes obey a set of coupled linear "Bogoliubov equations" 



Lu 3 - g (*) v 3 = Eju 



j > 



(57a) 



Lvj-g{**) 2 Uj=-EjVj, (57b) 

where 

L = T + V tT - f i + 2g\V\ 2 (58) 

is a Hermitian operator. 

Straightforward manipulations with the Bogoliubov equations show that Ej J dV {\uj\ 2 — |fj| 2 ) is real. If the 
integral f dV (\v,j\ 2 — |fj| 2 ) is nonzero, then Ej itself is real. Like Eq. ( [si]) for a uniform condensate, the Bose- 
Einstein commutation relations ( pa ) for the deviations from the nonuniform condensate can be shown to imply the 
following positive normalization |95[| 

"dV(\u j \ 2 -\v j \ 2 )=l. (59) 

For each solution uj , Vj with eigenvalue Ej and positive normalization, the Bogoliubov equations always have a 
second solution v*,u* with eigenvalue — Ej and negative normalization. The only exception to the requirement 
of real eigenvalues arises for zero-norm solutions with J dV (\uj\ 2 — \vj\ 2 ) = 0. In this case the character of the 
eigenvalue requires additional analysis. Numerical investigations |B7j of vortices in nonuniform trapped condensates 
have reported imaginary and/or complex eigenfrequencies for doubly quantized vortices but only real eigenfrequencies 
for singly quantized vortices. Specifically, for a repulsive interparticle interaction, Pu et al. fl97| found that singly 
quantized vortices are always intrinsically stable; in contrast, multiply quantized vortices have alternating stable and 
unstable regions with complex excitation energy as the interaction parameter Na/d increases. The most unstable 
vortex state decays after several periods of the harmonic trapping potential. In the case of multiply quantized vortices 
(q > 1), the vortex core contains localized quasiparticle bound states with small exponential tails; these modes have 
complex frequencies and are responsible for splitting the multicharged core |p8[ . For an attractive interaction, stable 
vortices exist only for the singly quantized case in the weak-interaction regime; a multiply quantized vortex state is 



always unstable. Similar imaginary and complex solutions have been found for dark solitons J99|-101|. For additional 



results on complex eigenfrequencies, see Ref. [102] and the Appendix of Ref. 1 103 1 
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In terms of the quasiparticle operators, the approximate perturbation Hamiltonian operator takes the simple intu- 
itive form 



K' ^Y^E j0 c\ aj , (60) 

j 

apart from a constant ground-state contribution of all the normal modes. Here, the sum is over all the states with 
positive normalization, and it is clear that the sign of the energy eigenvalues Ej is crucial for the stability. If one or 
more of the eigenvalues is negative, the Hamiltonian is no longer positive definite, and the system can lower its energy 
by creating quasiparticles in the unstable modes. 

The present derivation of the Bogoliubov equations and their properties emphasizes the quantum-mechanical basis 
for the positive normalization condition (^9|) and the sign of the eigenvalues. It is worth noting an alternative purely 
"classical" treatment [ 104| , 3~L| ] based directly on small perturbations of the time-dependent GP equation (^) around 



the static condensate ^(r). The solution is assumed to have the form 

*(r, t) = e~ iflt/n [*(r) + u{r)e~ luJt - v*(r)e lujt ] , (61) 
and the appropriate eigenvalue equations then reproduce Eqs. (Is?]). 



B. Uniform condensate 



For a uniform condensate, the solutions of Eq. (57) are plane waves, and the corresponding energy is the celebrated 
Bogoliubov spectrum |3q| 



E k = \J gntfk 2 /M + (h 2 k 2 /2M) 2 , (62) 

where k is the wave vector of the excitation and n is the condensate density. For long wavelengths fe(< 1, Eq. (|62| ) 
reduces to a linear phonon spectrum E k rs Hsk with the speed of compressional sound s = yj gn/M given by Eq. (|1C|). 
In the opposite limit k£ 3> 1, the spectrum reduces to the free-particle form plus a mean-field Hartree shift from the 
interaction with the background condensate E k ~ (H 2 k 2 /2M) + gn. 

To understand the importance of the sign of the eigenfrequency, it is instructive to consider the case of a condensate 
that moves uniformly with velocity vq. As noted in connection with Eq. (|24|), the condensate wave function is 
\£(r) = y / ne lq r , where q = Mvq/H and the chemical potential becomes fj, = \Mvq + gn. The Bogoliubov amplitudes 
for an excitation with wave vector k relative to the moving condensate have the form 



Uk(r) 
«k(r) 



(63) 



where the different signs ±iq • r arise from the different phases ±i2q • r in the off-diagonal coupling terms in the 
Bogoliubov equations (^)- The solution with positive norm has the eigenvalue 

E k (v ) = hk-v + E k , (64) 



as expected from general considerations [105,106]. In the long-wavelength limit, this excitation energy reduces to 
^k(vo) ~ hk(vocos0 + s), where 9 is the angle between k and Vo. For vq < s, the quasiparticle energy is positive 
for all angles 9, but for vq > s, the quasiparticle energy becomes negative for certain directions, indicating the onset 
of an instability. This behavior simply reflects the well-known Landau critical velocity for the onset of dissipation, 
associated with the emission of quasiparticles. It has many analogies with supersonic flow in classical compressible 



fluids [107] and Cherenkov radiation of photons in a dielectric medium [108,109]. For vq > s, the GP description 



becomes incomplete because the excitation of quasiparticles means that the noncondensate is no longer negligible. 

C. Quantum- hydrodynamic description of small-amplitude normal modes 

The quantum-hydrodynamic forms ( |l6| ) and ( |l7| ) of the time-dependent GP equation provide a convenient alter- 
native basis for studying the small-amplitude normal modes. The small perturbations in the density n'e~ lut and the 
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velocity potential $'e %UJt obey coupled linear equations |)6,110 1 1 1 1 that reduce to Eq. ( |33| ) in the TF limit for a 



static condensate | p7| . A comparison with Eqs. ( p6| ) shows that the quantum- hydrodynamic amplitudes 

n'j = ** Uj - * Vj = |*| {e~ iS Uj - e iS Vj) , (65a) 

$' = ($* u . + $ „ .) = f !— ( e - lS u, + e iS v~) (65b) 

are simply linear combinations of the Bogoliubov amplitudes Uj and Vj in the presence of the given condensate solution 
* = e l^l. The positive normalization condition (p9[) yields the equivalent quantum- hydrodynamic form 



JdVi(n^' j -&;n$ = ± (66) 
For many purposes, the quantum-hydrodynamic modes provide a clearer picture of the dynamical motion. 

D. Singly quantized vortex in axisymmetric trap 

Early numerical studies for small and medium values of the interaction parameter Na/do < 1 examined the small- 
amplitude excitations of a condensate with a singly quantized vortex |l!2| . In particular, the spectrum contained an 
"anomalous" mode with a negative excitation frequency and positive normalization associated with a large Bogoliubov 
amplitude u localized in the vortex core (see also relevant comments in Ref. p3| concerning the relationship between 
the sign of the normalization and the sign of the eigenfrequency) . The anomalous mode corresponds to a precession 
of the vortex line around z axis. As seen from the general discussion of the Bogoliubov equations, this anomalous 
mode indicates the presence of an instability. 

Since the condensate wave function has an explicit phase *(r) = e 1 ^ |*(rx, z)\, the Bogoliubov amplitudes for an 
excitation with angular momentum mh relative to the vortex condensate take the form 

u m (r) \ = f e^e^u m (r ± ,z)\ 
v m (r) J \e-^^v m (r x ,z)) ' 

analogous to those in Eq. ( |63| ) for a condensate in uniform motion. Here, the azimuthal quantum number m 
characterizes the associated density and velocity deformations of the vo rtex proportional to e m ^ [for example, 
n' m = 1*1 (u m — v m ) e lm ^ , as is clear from Eqs. (Bq)]. The numerical studies 112 found that the anomalous mode has 



an azimuthal quantum number m a — —I. Its frequency oj a is negative throughout the relevant range of Na/do < 1; 
in the noninteracting limit, u a approaches — u)j_, and cu a increases toward from below with increasing Na/do. 

To understand the particular value m a = —1, it is helpful to recall the noninteracting limit, when the negative 
anomalous mode for the vortex condensate signals the instability associated with Bose condensation in the first excited 
harmonic-oscillator state with excitation energy fku± and unit angular momentum. A particle in the condensate can 
make a transition from the vortex state back to the true harmonic-oscillator ground state, with a change in frequency 
— loj_ and a change in angular momentum quantum number —1. More generally, the density perturbation n' a for the 
anomalous mode with negative frequency — | is proportional to exp [i (\u a \t — <j))\ and hence precesses in the positive 



sense (namely counterclockwise) at the frequency \uj a 
the precession frequency of a one-component vortex [|3 



Thus the anomalous mode describes the JILA observations of 



1. Near-ideal regime 



An explicit perturbation analysis [113 68 of the GP equation for the condensate wave function in the weakly 



interacting limit found the thermodynamic critical angular velocity 



1 Na „ r ^,., ( Na^ 2 



wj_ V87T d z \ d z ' 



where the second-order correction depends explicitly on the axial asymmetry A = lu z /luj_. Similarly, a perturbation 
expansion of the Bogoliubov equations in the weak-coupling limit verified the numerical analysis and found the explicit 
expression for the frequency of the anomalous mode 
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-JL = -! + -—— + W W(A) -p +••• . (69) 
w_l V87r a z V <W 

It is evident that ft c + w Q vanishes through first order, and the detailed analysis shows that the second-order contri- 
bution to the sum is positive for all values of the axial asymmetry parameter A. 

The physics of the anomalous mode can be clarified by considering an axisymmetric condensate in rotational 
equilibrium at an angular velocity ft around the z axis. In the rotating frame, the Hamiltonian becomes H — QL Z , 
and the Bogoliubov amplitudes have frequencies Wj(ft) = ujj — rrijil, where u>j is the frequency in the nonrotating 



frame and rrij is the azimuthal quantum number [see Eq. (67)]. For the anomalous mode with m a = — 1, the resulting 
frequency in the rotating frame is 

W o (0) = 0J a + ft, (70) 

which is directly analogous to Eq. ( |64| ) for uniform translation. Since ui a is negative, the anomalous frequency in 
a rotating frame increases linearly toward zero with increasing ft; in particular, u) a (Q) vanishes at a characteristic 
rotation frequency 

ft* = -cj a = K| (71) 

that signals the onset of the regime ft > ft* for which the singly quantized vortex becomes locally stable. Equation ( |69| ) 
gives an explicit expression for ft* in the weak-coupling limit, and detailed comparison with Eq. (^8|) indicates that 
ft* < ft c for any axial asymmetry A (but only because of the second-order contributions). It is natural to identify 
ft* with the angular velocity for the onset of local stability with respect to small perturbations; this quantity was 
denoted ft m in connection with the equilibrium energy in the TF limit (see Fig. ||). 



2. Thomas-Fermi regime for disk-shape trap 



The anomalous negative-frequency mode exists only because the condensate contains a vortex. Hence it cannot be 
analyzed by treating the vortex itself as a perturbation. In the TF limit, however, it is possible to use Gross's and 
Pitaevskii's ||Q solution (|l|) for a vortex in a laterally unbounded fluid as the basis for a perturbation expansion. 
A detailed analysis of the Bogoliubov equations for an axisymmetric rotating flattened trap in the TF limit yields the 



explicit expression for the anomalous mode [114 



As in Eq. (fry) for the weak-coupling limit, Eq. d72) yields 



(72) 



ft* = i 



2 MRj 



■In 



R± 



5 



(73) 



where the last two equalities follow from ( |5l| ) and ([52]). This relation further supports the identification of ft* with 
the metastable rotation frequency ft m associated with local stability of a vortex for small lateral displacements from 
the center of the trap. Note that ft m < ft c for a disk-shape condensate (in the TF limit) [see Eqs. (|5l| ) and (M)], 
similar to the behavior for the weak-coupling regime. 



3. Quantum-hydrodynamic analysis of condensate normal modes in the Thomas-Fermi regime 

In addition to the anomalous mode described above, the condensate has a sequence of normal modes that occur 
both with and without a vortex. Indeed, one of the early triumphs of the quantum-hydrodynamic description ]57| was 
the detailed agreement between the theoretical predictions and the measured frequencies of the lowest few collective 
normal modes pl[ . For an axisymmetric condensate, the normal modes can be classified by their azimuthal quantum 
number m, and modes with ±m are degenerate for a stationary condensate. 

When the condensate contains a vortex, however, the various collective modes are perturbed. In particular, the 
vortex breaks time-reversal symmetry by imposing a preferred sense of rotation, so that modes with ±to are split 
(this behavior is analogous to the Zeeman effect in which an applied magnetic field splits the magnetic sublevels) . In 
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FIG. 11. Transverse oscillations of a stirred condensate with 3.7 x 10 5 atoms in an elongated trap with u±/2ir = 171 Hz. 
For (a)-(c), the stirring frequency Q,/2n = 114 Hz is below the threshold for vortex nucleation, whereas for (d)-(f), the stirring 
frequency Q/2n — 120 Hz has nucleated a vortex (visible at the center of the condensate). The sequences of pictures correspond 
to time delays r = 1, 3 and 5 ms for which the ellipticity in the xy plane is maximum. The fixed axes indicate the excitation 
basis of the quadrupole mode and the rotating ones indicate the condensate axes. 
(Taken from Ref. f§). 



fact, the splitting of these degenerate hydrodynamic modes has been used to detect the presence of a vortex |B6 92 
and to infer its circulation and angular momentum. 

In the context of the quantum-hydrodynamic description, the principal effect of the vortex arises through its 
circulating velocity field v, which shifts the time derivative dt — ► dt + V • v. For a normal mode oc e %m ^ with 
azimuthal quantum number m, the perturbation in the frequency has the form mhjMr\. A detailed analysis shows 
that the fractional splitting of the modes is of order (uj + — u)-)/lo + ~ |m|dj_/i?j_, with a numerical coefficient that 
depends on the particular mode in question f62,lll|. Independently, Zambelli and Stringari 115 used sum rules to 
calculate the vortex-induced splitting of the lowest quadrupole mode with m = ±2; the two approaches yield precisely 
the same expressions. In the absence of a vortex, the \m\ = 2 mode simply involves an oscillating quadrupole distortion, 
but the vortex-induced splitting means that the quadrupole distortion precesses slowly in a sense determined by the 
circulation around the vortex. The angular frequency of precession of the eigenaxes of the quadrupole mode is equal 
to (u> + — u)J)/2\m\ = (uj+ — w_)/4 = Xoj±dj_/Rj_. Figure |ll| shows the difference between the two cases (with and 
without a vortex) for a condensate with w 3.7 x 10 5 87 Rb atoms in an elongated trap with ui±/2ir = 171 Hz. In 
the ENS experiment |36j ], when one vortex is nucleated at the center of the condensate, the measured frequency 
splitting of the quadrupole mode (lu + /2tt — 250 Hz) is (w+ — oj_)/27t = 66(±7) Hz. For the experimental parameters 
(R± = 3.8 /urn), theory predicts (uo + — lu-)/2ii — 7h/2nMR\ = 56 Hz. The result holds in the TF limit and is valid 
with an accuracy of order d\ ln(R± J '£) / 'R\ ~ 0.15. With this uncertainty, the theoretical prediction 56(±8) Hz agrees 
with the experimental value. 

One should note that the vortex-induced splitting of the condensate modes is maximum if the vortex is located at 
the trap center. If a straight vortex line is displaced a distance £o = tq/R± from the z axis of the TF condensate, 
then the splitting of the quadrupole mode (m = ±2) is given by the expression 



ijJ+ — u>- = lu>\ 



1 - ~fi 



i 4- -r 4 - r 6 + — r 8 



(74) 



The splitting goes to zero if the vortex moves out of the condensate (Co = !)• 
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4- Numerical analysis for general interaction parameter 



Garcfa-Ripoll and Perez-Garci'a |102|| have performed extensive numerical analyses of the stability of vortices in 
axisymmetric traps with an axial asymmetry parameter A = oj z /oj± — 1 (a sphere) and A = i (one particular cigar- 
shape condensate). They conclude that a doubly quantized vortex line has normal modes with imaginary frequencies 
and that an external rotation cannot stabilize it. For a singly quantized vortex in a spherical trap, however, they 
confirm the presence of one negative-frequency (anomalous) mode with \ui a \ < 51 c . For their cigar-shape condensate, 
they find additional negative-frequency modes and suggest that such elongated condensates are less stable than 
spherical or disk-shape ones. More recent numerical work 116 83] confirms these findings for other geometries, 
especially that for the ENS experiment (34), where the axial asymmetry is large (u>±/u; z ~ R z /R± ~ 14). It is 
expected that a vortex in an elongated condensate becomes stable only for an external angular velocity Q m = max |ct> |> 
where max|w | is the absolute value of the most negative of these anomalous modes. For only modestly elongated 
traps, the metastable frequency fl m exceeds the thermodynamic critical value fi c ; these results provide an alternative 
explanation of the ENS observation that the first vortex appears at an applied rotation ss 70% higher than fl c . 
Independently, an analysis of the bending modes of a trapped vortex Jll7[ in the TF limit finds that a vortex in a 
spherical or disk-shape condensate has only one negative frequency (anomalous) mode, but the number of such modes 
in an elongated condensate increases with the axial asymmetry ratio R z /R± (discussed below in Sec. V.D.2). 



V. VORTEX DYNAMICS 



The preceding sections considered the equilibrium and stability of a vortex in a trapped Bose condensate, using 
the stationary GP equation and the Bogoliubov equations that characterize the small perturbations of the stationary 
vortex. These approaches are somewhat indirect, for they do not consider the dynamical motion of the vortex core. 
The present section treats two different methods that address such questions directly. 



A. Time-dependent Variational Analysis 



Consider a variational problem for the action J dt C(t) obtained from the Lagrangian 



C{t) 



dV 



2 I dt dt 



**(T + v tI -nL z )y-y\*\ 



(75) 



It is easy to verify that the Euler-Lagrange equation for this action is precisely the time-dependent GP equation in 
the rotating frame. 

If, instead of ^t(r,t), we substitute a trial function that contains different variational parameters (for example, 
the location of the vortex core), the resulting time evolution of these parameters characterizes the dynamics of the 
condensate. This method is not exact, but it provides an appealing physical picture . For ex ample, it determined 
the low-energy excitations of a trapped vortex- free condensate at zero temperature | 11S| , 11S| for general values of 
the interaction parameter. In the TF limit, this work reproduced the expressions derived by Stringari |57|] based on 
Eq. (1). 



1. Near-ideal regime 

In the near-ideal limit, only the axisymmetric case has been studied, and it is natural to start from the noninteracting 
vortex state (^6|), incorporating small lateral displacements of the vortex and the center of mass of the condensate, 



along with a phase that characterizes the velocity field induced by the motion of the condensate [ 12C ] . In addition to 
the rigid dipole mode (in which the condensate and the vortex oscillate together at the transverse trap frequency 
an extra normal mode arises at the anomalous (negative) frequency ui a given in Eq. (|6^) omitting the second-order 
corrections that are beyond the present approximation. In this weak-coupling limit, the resulting displacement of the 
vortex is twice that of the center of mass, so that both must be included to obtain the correct dynamical motion. 
Detailed analysis confirms the positive normalization and relative displacements found from the Bogoliubov equations 



for the same axisymmetric trap [113] 
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2. Thomas-Fermi regime for straight vortex in disk-shape trap 



For a nonaxisymmetric trap in the TF regime, only the nonrotating case (f2 = 0) has been analyzed, using the fully 
anisotropic TF wave function as an appropriate trial state, again with parameters describing the small displacements 
of the straight vortex and the center of mass of the condensate ]7^] . The trial wave function was chosen in the form 

* (r,i) = B(t)f [r - r (t)] F[r - T] (t)] ]J exp [ix 3 a 3 {t) + ix]f3 3 {t)] . (76) 

j=x,y,z 

Here the function / (r) characterizes the vortex line inside the trap and far away from the vortex core has the 
approximate form / (r) = e l ^\ the function F(r) is the TF condensate density. The time-dependent vector r) (t) = 
(iJoxiVoy-iVoz) describes the motion of the center of the condensate, while ro{t) = {xo,yo,0) describes the motion 
of the vortex line in the xy plane. The other variational parameters are the amplitude B(t) and the set aij(t) and 
Pj{t). Substitution of the trial wave function into ((75|) yields an effective Lagrangian as a function of the variational 
parameters (and their first time derivatives). The resulting Lagrangian equations have a solution that corresponds to 
the motion of the vortex relative to the condensate. For this solution the vortex motion is described by 

x a = s Q R x sin (u a t + (f>o) , y Q = e R y cos (uj a t + cf) Q ) , (77) 

while the displacement of the condensate is given by 

15e £ 2 fR±\ R x 
Vox = 2R^ \T ) R~TR~ Sm 



15e o£ 2 ,_{R±\ Ry 



where 



m y = — ln \~T ) Rx + Ry CQS ^ at + ^ ' ^ 



3fu, x u y ln / R±\ = 3H ln /J^ (8Q) 



4/x V £ / 2MR x R y V C 



in agreement with that found in Eq. (72). The quantity Xq jR x -I- j/q / Ry — e$ remains constant as the vortex line follows 



an elliptic trajectory around the center of a trap along the line Vt r = const, and the energy of the system is conserved 



[as follows from Eq. (49)]. The condensate also precesses with the relative phase shift ir at the same frequency, but 



the amplitude of the condensate motion is smaller than that of the vortex line by a factor ~ £ In (Rj_/\q\£) /R X R 



For an axisymmetric TF condensate in rotational equilibrium at an angular velocity f2, the Lagrangian ( |T5|) provides 
a more general result for the precession frequency. With the hydrodynamic variables "J = e*' s |^|, the first term of 
the Lagrangian becomes —H J dV \"i>\ 2 dS/dt. Since the TF condensate density vanishes at the surface, the particle 
current also vanishes there, and it usually suffices to assume a single straight vortex displaced laterally to ro(t), with 
S(r, r ) = arctan[(y — yo)/(x — xq)] and no image vortex. The Lagrangian becomes 

£ = J dVMn(r)r -Mr)-E(r ) + nL z (r ), (81) 



where 



v (r) = A V S(r,r ) = -Av S(r,r ) = ( K /2n) ^ g> (82) 
M M r — rp 



is the circulating velocity field about the vortex line. In the special case of a two-dimensional condensate with the 
TF density n(r) — n(0)(l — r\/R 2 ^) per unit length, Eq. (|8ll ) becomes 

C 2 = {h + tyL z2 (r ) - <p Q L z2 (0) - E 2 (r ), (83) 
where <fro is the azimuth angle describing position of the vortex line, 

L Z 2(r ) = \n{0)nRlh(l - C 2 ) 2 (84) 
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and 



E 2 (r ) 



K 2 Mn(0) 



2(1 -^.ln(^)+(1-C,T)lu(:i -„T)- 1 + 2C„ 



(85) 



with Co — 1*0 /R± [note that ^(0) is the mean particle density n per unit length]. These expressions differ from the 



classical results for a uniform fluid in a rotating cylinder |88 121] because of the parabolic TF density; in particular, 
the TF angular momentum per unit length L z2 here is proportional to (1 — Co) 2 j whereas that for a uniform density 
is proportional to 1 — £q. 

The Lagrangian dynamical equations show that the vortex precesses at fixed tq with the angular frequency 

^ = -n+ dE2/dr ° = _n - dE2/dr ° (86) 

dL z2 /dr KMr n(r ) ' 



This result is just that expected from the Magnus force on a straight vortex [ 122j - 124] . For small displacements from the 



center, the precession frequency in a nonrotating two-dimensional condensate reduces to 4>q ~ (k/2ttR 2 _) ln(i?j_/£) 
i r2 c | ]70| , but </>o increases with increasing r and eventually diverges near the edge of the condensate where the density 
vanishes. 

The corresponding results for a three-dimensional disk-shape TF condensate follow from Eqs. ( fl9| ) and (pl|). In 
particular, the integration over z means that the total angular momentum L z % = Nh(l — Co) 5 ^ 2 associated with the 
presence of the vortex differs from the two-dimensional result proportional to (1 — Co) 2 - Apart from numerical factors 
reflecting the three-dimensional geometry, Eq. (^) remains correct. For a straight vortex, it yields 

where fl m = ^(H/MR t j_)hi(R±/{;) is the metastable frequency ( |5~ij ) for the appearance of a central vortex in a disk- 
shape condensate. In the special case of a vortex near the center (tq — > 0), this precession frequency reduces to 
(minus) the corresponding anomalous frequency w a (f2) in Eq. ([72]) for a condensate with a single central vortex line. 
To understand why the precession frequency <^o is the negative of the anomalous frequency, recall that the linearized 
perturbation in the density for the anomalous mode is proportional to expi[m a — u> a (£l)t] = expi[— <f> — u> a (Cl)t] 
because m a = — 1; this latter form shows clearly that the normal mode propagates around the symmetry axis at an 
angular frequency — u> a (Q), with the sense of rotation fixed by the sign of — Lu a (Q). 

According to (^7|), for a nonrotating trap the precession velocity of a displaced vortex increases with the vortex 
displacement as v = O m ro/(1 — r{\/R\). It is interesting to estimate at what displacement the vortex velocity becomes 
supersonic |l25| . Assuming the speed of sound varies radially with the local density as c = Co\fl — r\jR\, where 
c = \J Ji/M = luj_R±/\/2, we obtain v/c = (\/2^mTo/^±-R±)(l — Tq/ 'R\)~ 3 / 2 . As a result, the vortex velocity 
becomes supersonic if 

" >^U-4r) m = ^^('-4r) m . (83) 



R± V2n m V R±J 3£ln(i? x /£) V R 



For parameters of JILA experiments |37| i?j_/C" « 33, this gives a critical displacement of r^/R^ sa 0.82 where the 
precession vortex velocity becomes supersonic. 



B. Method of Matched Asymptotic Expansions 

At zero temperature, the dynamics of a condensate in a rotating nonaxisymmetric trap follows from the appropriate 
time-dependent GP equation 

d*t / S 2 V 2 o \ 

m ~m = \ W + VtI + m ~ ll{sl) + iMl ■ (r x V) J *' (89) 

A vortex line in the condensate will, in general, move in response to the effect of the nonuniform trap potential and 
the external rotation, as well as self-induced effects caused by its own local curvature. This problem can be solved 
in the case of a large condensate, where the TF separation of length scales means that the vortex-core radius £ is 
much smaller than the condensate radii Rj. The relevant mathematics involves the method of matched asymptotic 



expansions [ 126 - 128 1 
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1. Dynamics of straight vortex in Thomas- Fermi regime for disk-shape trap 



As an introduction to these techniques, it is helpful first to concentrate on the case of a straight singly quantized 
vortex line Q , which is applicable to disk- shap e condensates with R z <C R± ; this analysis generalizes two-dimensional 
results found by Rubinstein and Pismen ]l27f| . Assume that the vortex is located near the center of the trap at a 
transverse position r±o(t). In this region, the trap potential does not change significantly on a length scale comparable 
with the vortex core size £. The method of matched asymptotic expansions compares the solution of Eq. ( |S9| ) on two 
very different length scales: 

First, consider the detailed structure of the vortex core. Assume that the vortex moves with a transverse velocity 
Viz, and transform to a co-moving frame centered at the vortex core. Away from the trap center, the trap 
potential exerts a force proportional to Vj_T4 r evaluated at the position r±o(t). The resulting steady solution includes 
the "asymptotic" region |rj_ — rj_ | S> £. 

Second, consider the region far from the vortex (on this scale, the vortex core is effectively a singularity). The 
short-distance behavior of this latter solution also includes the region £ -C |rj_ — rj_o|. The requirement that the two 
solutions match in the overlapping region of validity determines the translational velocity V of the vortex line. 

Unfortunately, the details become rather intricate, but the final answer is elegant and physical: 



v= m 



4Mfx 



hi 



r ± \ 



(z x V_LV tr ) = 



(z x V ± V te ) , (90) 



where Rj_ for an asymmetric trap is defined in Eq. (50). This expression has several notable features. 

(a) The motion is along the direction z x Vj_Vt r and hence follows an equipotential line of Vt r . Thus the trajectory 
conserves energy, which is expected because the GP equation omits dissipative processes. In the present case of an 
anisotropic harmonic trap, the trajectory is elliptical. 

(b) For a nonrotating trap (Q — 0), the motion is counterclockwise in the positive sense at the frequency given by 
Eq. (pTj|), proportional to io x u y . 

(c) With increasing applied rotation f2, the translational velocity V decreases and vanishes at the special value 



^-^t-^if r^nf (91) 

proportional to | (ur%. + ujy) . This value precisely reproduces Eq. (|5l]) associated with the onset of metastability for 
small transverse displacements of the vortex from the trap center. 

(d) For £1 > fi OT) the motion is clockwise as seen in the rotating frame. A detailed analysis based on the normalization 
of the Bogoliubov amplitudes shows that the positive-norm state has a frequency [compare Eq. (|8fj|)l 

Wa (n ) = 4^L(0-fi ro ). (92) 

Note that this expression differs somewhat from Eq. ( |70| ) because the trap here is anisotropic. The normal-mode 
frequency is negative and hence unstable for Q < J7 m , but it becomes positive and hence stable for Q > VL m . 

This direct analysis of the motion of a straight vortex reproduces the physics of the onset of (static) metastability 
( pf| ) studied with the GP Hamiltonian and the (dynamic) anomalous mode ( f73| ) and (]8(]) studied with the Bogoliubov 
equations and with the Lagrangian method. 



2. Dynamics of curved vortex in Thomas-Fermi regime 

Consider a nonaxisymmetric trap that rotates with an angular velocity H (for convenience, f2 is often taken along 
the z axis). At low temperature in a frame rotating with the same angular velocity, the trap potential is time 
independent, and Eq. (^9|) describes the evolution of the condensate wave function. In the TF limit, the method of 
matched asymptotic expansions again yields an approximate solution for the motion of a singly quantized vortex line 
with instantaneous configuration ro(z, t). Let t be the local tangent to the vortex (defined with the usual right-hand 
rule), h be the corresponding normal, and b = t x n be the binormal. A generalization of the work of Pismen and 
Rubinstein [126,127] eventually yields the explicit expression for the local translational velocity of the vortex [117] 



27 



h f ixvv t ,(r a ) a / i , fc^ , 2vy tr (r )xn 



where fe is the local curvature (assumed small, with /;( C 1) and Aj_ is the Laplacian operator in the plane perpen- 
dicular to SI. 

This vector expression holds for general orientation of the gradient of the trap potential, the normal to the vortex 
line, and the angular velocity vector. Near the TF boundary of the condensate, the denominator of the first term 
becomes small, implying that the numerator t x V Vtr(ro) must also vanish near the boundary. As a result, the axis 
of the vortex line t is parallel to VVtr at the surface and hence obeys the intuitive boundary condition that the vortex 
must be perpendicular to the condensate surface. 



C. Normal modes of a vortex in a rotating two-dimensional TF condensate 



This very general Eq. ( |93| ) applies in many different situations [117]. The simplest case is an initially straight vortex 
in a two-dimensional asymmetric TF condensate with SI = Slz and u> z = (hence no confinement in the z direction). 
For small displacements, the x and y coordinates of the vortex core execute harmonic motion oc exp[i(nz — tot)] that 
can vary between helical and planar depending on the relative phase of the x and y motion. The dispersion relation 
lu k (SI) depends on the continuous parameter k and the rotation frequency SI, along with the TF radii R x and R y 



where 



"'" (n) = ± uhr, V ( 2 " - ") ( 2 - " 2fl » " & ) ta War + « ' (94) 



AMR 2 x Rl 

SI = x y —rSl (95) 



H(Rl + Rl)ln^^ + if 

is a dimensionless rotation speed. 

Of all the various normal modes, a straight vortex line (k = 0) has the most negative (anomalous) frequency 



w„(0) = - 



2MR x Ry 



ln 



R± \ 4^Sl 



(96) 



where an analysis similar to that for Eq. (92) shows that the minus sign corresponds to the Bogoliubov solution 
with positive norm. For O = 0, the vortex precesses counterclockwise about the z axis in the positive sense. With 
increasing rotation frequency SI, the precession frequency decreases and vanishes at SI = Sl m , where the metastablc 
rotation frequency in two dimensions is 

h{ul + ul) fR ± \ h (Rj_\ 

as expected, this value is the precession frequency ^Sl c discussed below Eq. (^6|) [compare Eq. ( pl| ) for Sl m in a 
three-dimensional disk-shape TF condensate; the different numerical coefficient arises from the integration over the 
parabolic density in the z direction]. 

More generally, for n 2 > and a nonaxisymmetric trap (R x > R y ), the oscillation frequency can be imaginary 



(and hence unstable) within a range of axial wave numbers determined by y (2 — Sl)/R x < \n\ < y (2 — St)/Ry. For 
sufficiently fast rotation, however, the frequencies become real, and the small oscillations become stable at a rotation 
frequency SI > Sl m = 2. In the limit of a uniform unbounded condensate (R x ,Rj, — > oo), the general dispersion 
relation reduces to the familiar one for helical waves on a long straight vortex line 

c = ±A K 2i n(N0 . (98 ) 
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Using this dispersion relation, Barenghi |129[ estimated the amplitude of the vortex waves due to thermal excitation 
(the cloud is assumed to rotate at an angular velocity £1 > 57 c , so that the vortex is stable). He showed that finite- 
temperature effects in a Bose condensate can distort the vortex state significantly, even at the very low temperatures 
relevant to the experiments. For T = K, n w 10 12 — 10 13 cnT 3 and R w 5 /im, the amplitude of vortex oscillations 
can be 4-14 times the size of the vortex core. At the same time, the thermal excitation of vortex waves in superfluid 
4 Hc is negligible (much smaller than the corresponding vortex-core size). 

D. Normal modes of a vortex in a rotating three-dimensional TF condensate 

Consider a three-dimensional TF condensate with lo z > 0, confined within a TF region z 2 < R 2 = 2/z/Afw 2 . 



1. General formalism 



For a vortex that initially lies along the z axis, it is straightforward to find the pair of coupled equations for the 
small transverse displacements of the vortex x(z,t) and y(z,t). In particular, we seek solutions of the form 



x — x(z) sin(tjt + ipo), y = y(z) cos(wi + ipo), 



(99) 



in which case the amplitudes x(z) and y{z) describe the vortex shape and obey coupled ordinary differential equations. 
Introducing dimensionless scaled coordinates x — > R x x, y — > R y y, z — > R z z, we find from Eq. (93) 



w(l — z 2 )x = — — 
dz 



y + n(l-z 2 )y, 



(100) 



cD(l 



z 2 )y 



where 



dz 



a{l-z 2 ) 



dx 
dz 



+ £l(l -z 2 )x, 



(101) 



(102) 



characterize the trap anisotropy and 



2MR x R y 
n \n(R±/0 



uj, Q 



4MR 2 x R 2 y 



h(R 2 x + R 2 y )HR ± /0 



(103) 



are dimensionless angular velocities. 



These equations (100) and (101) constitute a two-component Sturm-Liouville system with natural boundary con- 



ditions |130) | because the factor 1 — z vanishes at z — ±1. Consequently, the eigenfunctions merely must remain 
bounded at the surface of the condensate. A straightforward generalization of the usual analysis shows that the 
eigenfunctions obey the orthogonality condition 



dz (1 - z 2 )x m y n oc 8 n 



(104) 



2. Special solutions 

In the general case of a nonaxisymmetric trap, the resulting equations remain coupled, but they separate in the 
particular case of stationary solutions with u) = 0. For a nonrotating trap, such configurations reflect a balance 
between the effects of curvature and the nonuniform trap potential. For example, the small-amplitude stationary 
solutions x n {z) remain finite at the surface z = ±1 only for certain special values of the trap anisotropy 
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where n > is an integer. The corresponding solutions have the form x n (z) oc P n (z), where P n is the familiar Legcndre 
polynomial. The solutions have n nodes and cross the z axis n times. If a differs from one of these special values 



(105), there is no stationary configuration. Similarly the equation for the y displacement has stationary solutions 
only if (3 = R 2 y /R 2 z = 2/[m(m + 1)]. 

This classification of the solutions by the number of nodes remains more generally valid. In the special case of an 
axisymmetric condensate (a = /?) , we can consider the precession frequency ui n of the mode with n nodes as a function 
of the axial trap anisotropy a. Evidently the function iv n changes sign at the special value a = a n = 2/[n(n+l)]. This 
observation allows us to determine the number of modes with negative frequencies at a fixed value of the anisotropy 
parameter a. For a > 1 (a spherical or disk-shape condensate), only one mode has a negative frequency. If ^ < a < 1, 
there are two such anomalous modes, and so on. If a n < a < a nonrotating axisymmetric TF condensate has 

n anomalous modes with negative frequency. 

The special case of a nearly disk-shape anisotropic rotating TF condensate is particularly tractable because a -1 
and (3~ x provide small expansion parameters. There is only one relevant normal mode, with frequency 

w„(fi) = -fi m + 0, (106) 



where 



V I 



8// 



— (- - 

10 \a + [3 



In ( ) lor nearly disk-shape TF condensate. (K)i ) 



If Q < Q m = |w a (0)Lthe frequency is negative, and the mode is therefore unstable. This value generalizes that found 
previously in Eqs. (|5l| ) and ( |9l| ) for the angular velocity at which a straight vortex at the center of a thin disk-shape 
condensate becomes metastable, now including the first corrections of order cT 1 and /3 _1 . 

This result ( |107|) remains approximately correct for a spherical condensate [a, = j3 = 1), which is the geometry 
used in recent JILA experiments |[37fl. Since f2 m is numerically equal to the frequency \ui a \ of the one anomalous 



mode in the nonrotating condensate, Eq. (107) also yields the precession frequency of a nearly straight vortex moving 
counterclockwise around the center of the condensate [ |78|]83| . In particular, we find \w a \/(jj — f(£/i?) ln(1.96i?/£), 
where u> is the isotropic trap frequency and the additional numerical factor 1.96 in the logarithm is the next correction 
to the logarithmic accuracy (see, for example, Ref. |127| ]). With the JILA parameters R » 22 /im and £ » 0.67 /xm, 
this expression yields |w |/27r = 1.58 ± 0.16 Hz, where the uncertainty reflects the omission of corrections of relative 
order (£/R) ln(i?/£) w 0.1. For comparison, the experimental value 1.8 ± 0.1 Hz for the precession frequency is 
somewhat larger, but the theoretical prediction is sensitive to the number N of atoms in the condensate and, as seen 



in Eq. (BTl), to the radial displacement of the vortex [183,123 



The situation is very different for an elongated cigar-shape condensate with R z ^> R±, when the solutions for 
the precessing normal-mode amplitudes grow exponentially with \z\. In contrast to the two-dimensional case, such 
solutions are now possible because the condensate is bounded along the z axis. In the simplest case of an axisymmetric 
trap with R x = R y = R±, the mode with no nodes has a frequency w a (f2) = — fi m + f2. Although this expression has 



the same form as Eq. (106) for a disk-shape condensate, the physical behavior is very different because the metastable 
angular velocity 

n m = — ^5- M-]n( —) w fl c (108) 
2MR\ R\ V W 5R 2 ± v 1 

beco mes large for a highly elongated TF condensate. For the ENS geometry Q,^6|, where uj±/uj z R z /R± w 14.4, 



Eq. (lQg) is far too large to fit the observations and can even exceed the limit of rotational mechanical stability 



Q = uj±_ that occurs when the centrifugal force cancels the confining trap potential. 

For a harmonic transverse external potential oc r\, the method of matched asymptotic expansions is valid if the 
vortex displacement r from the z axis satisfies the condition r > £ (in the vicinity of the vortex core the trap potential 
is approximated as a linear function) . For a long cigar-shape condensate, the solution for the lowest mode has the form: 
r = r Q cosh(z/a), where r is the vortex displacement at z = 0. The condition of small vortex displacement implies 
that r cosh(l/a) -C R±, while the condition of small vortex curvature k£ <C 1 implies that r £ cosh(l/a)/i?2 a2 ^ 1- 



A combination of these conditions gives the following restriction on the validity Eq. (108): exp(l/a) <C 2R±/£. For 
the ENS experiments, 1/a ~ 200 and ~ 21, so that this condition fails. 
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1 2 3 4 5 

FIG. 12. Dimensionless frequencies u) = ui(Q = 0) for the first three normal modes of a vortex in an axisymmetric trap as 
a function of the axial anisotropy a — R^_/R%. The lower horizontal line is the negative of the dimensionless thermodynamic 
critical angular velocity Q c — 5. Note that |tDo| > for a < 0.26. 



As mentioned in Sec. IV. D. 4, the frequency for the onset of metastability ft m in Eq. (108) can be larger than 
the thermodynamic critical angular velocity ft c in Eq. (fl3|). This behavior is readily understandable because ft c 
characterizes the energy of a straight vortex along the symmetry axis [compare Eq. (Ejh] , whereas the most unstable 
normal-mode amplitude explicitly involves the small-amplitude distortion with no nodes. For a very elongated con- 
densate, the resulting vortex dynamics is particularly sensitive to the large curvature of the condensate surface near 
the two ends of the symmetry axis (in contrast to the small curvature for the flattened condensate) . 

Recent numerical studies [116,33| of the most negative anomalous modes for a trap geometry corresponding to 
the ENS experiments |54|,[56| yield values of ft m that are significantly smaller than the prediction given in Eq. (108). 
Reference ]116[] mentions the possible failure of the TF picture in the transverse direction, even though the conventional 



TF ratio is large, at least near the plane z = 0. As confirmation of the validity of the GP equation and the 

particular role of the anomalous modes, the numerically determined (8^] ft TO /27r » 0.73v± w 124 Hz agrees well with 
the ENS value ft b s /27r w 120 Hz for the appearance of the first vortex. 

For an axisymmetric trap (a = /?), we can seek normal-mode solutions in the form x(z) = y(z), leaving a single 
equation 



w(ft) - ft (1 - z 2 )x = - 



d 

dz 



a(l - z 2 ) 



dx 
dz 



(109) 



that depends only on the Doppler-shifted frequency w(ft) — ft = Q(0). The eigenfunctions are even or odd functions 
of z and can be classified by the number of times the vortex crosses the z axis (the number of nodes), m = 0, 1, 2, • • •. 
Figure 12 shows the dimensionless frequency u)(0) as a function of the trap anisotropy a = R\/R 2 for m = 0, 1, and 
2. In agreement with the analytical results, a disk-shape trap (a > 1) has only a single mode with negative frequency 
u>q. For -| < a < 1, there are two such modes (m — and m = 1) and successively more negative- frequency modes 
appear for smaller a. As noted previously, the critical frequency ft m for metastability is |£>o|, which is smaller than 
ft c for disk-shape traps and for moderately elongated traps. Our numerical analysis for the present TF limit predicts 
that ft m > ft c for a = R\/ R 2 Z < 0.26, which is somewhat larger than the value 0.2 implied by the limiting expression 



in Eq. (|08|). 

As in the case of a two-dimensional condensate, the frequency of the anomalous modes can become imaginary for 
an anisotropic trap with R x ^ R y [117]. To demonstrate that result, let us consider Eqs. ( |100[ ) and (101) for a trap 
close to axisymmetric with \a — f3\ -C a. The eigenfrequencies of the axisymmetric trap [with a = (3 = ao = h(a + (3)] 
are real and have the form <D TO (ft) = uj rn + ft, where m = 0, 1, 2, • • • denotes the various modes. For an anomalous 
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mode, the frequency Q m is negative, and the eigenfrequency oj m ((l) is equal to zero if the trap rotates with the angular 
velocity (l = \ui m \. One can rewrite Eqs. (100) and (101) as follows: 



where 



Cj{\-z 2 ) 



Ho 



V 



H = {-2 - a d z [(l - z 2 )8 z ] + (1 - z 2 )\u m \) f J 



1 




(110) 



V = -d z [(l-z 2 )d z 





a — ao 



/3 - a 




+ (i-z 2 )(n -\u n 



1 

1 



Considering V as a perturbation, we obtain the following expression for the normal-mode frequency in a nonaxisym- 
metric trap: 



u> = ± 



n 



where 



I-idz(l - z 2 ) (d z x m f 

J-rn — i > 

2/^(1- 



(111) 



(112) 



and x m — x m (z) describes shape of the mth vortex mode. As we increase the trap rotation, the eigenfrequency is 
real for ft < \uj m \ — I m \a — j3\ . Then, when |fi — \uj m \\ < I m \a — j3\, the frequency becomes imaginary. Finally, if 
ft > \u) m \ + I m \a — j3\, the frequency again becomes real. For a given trap anisotropy (given a and /3), one or several 
normal modes of the vortex have negative frequency. Trap rotation Q shifts the frequencies in the positive direction. 
When the frequency of a normal mode in the rotating frame approaches zero, the frequency becomes imaginary until 
\Q m = I m \a — f3\. If we increase the trap rotation further, the frequency (in the rotating frame) becomes positive. 
For a disc-shape condensate (with ao S> 1) there is only one anomalous mode with x a = y a = e (l + z 2 /2a$) and 



= -3 



. For a nonaxisymmetric rotating trap, the frequency of this mode becomes imaginary in the interval 



\R X — R y \/R x is the trap anisotropy in the transverse direction. Thus for a 
the solution has an imaginary frequency in a relatively narrow range of trap 



|fi - \Q a \\ < e|a) |/(15ao), where e = 
disk-shape condensate (with ao > 1 
rotation. 

For a cigar-shape condensate, several normal modes have negative frequencies. In the limit ao -C 1, the solution for 
the lowest anomalous mode has the form x a = y a = £cosh(2:/ao) and u> a ~ — l/ao- Consequently, the frequency is 
imaginary if |f2 — |o) a || < e|a) a |, namely in a relatively wide range of trap rotation. If the transverse trap anisotropy is 
large enough, several different anomalous normal modes can have imaginary frequencies in the same range of angular 
velocities. In this case a vortex along the z axis is stable (there are no normal modes with imaginary frequencies) only 
if the trap rotates slightly faster than the frequency of the lowest anomalous mode. This behavior could be relevant 
to ENS experiments. 



3. Energy of a curved trapped vortex 

Consider a trap that contains a singly quantized vortex and rotates with angular velocity f2 about the z axis. At 
zero temperature, Eq. ( |93"|) governs the dynamics of each element of the line 



_ r , , H AxWtr(r) , f V /, 1 k 2 \ 2W tr (r)xn , oS 

V < r > = -^bi^ + t T"^VV^) + A.vjr) ■ < 113 > 

where r = (x(z), y(z), z) determines the shape of the line. Correspondingly, Eq. (flCj) serves as the energy functional 



= JdV (||j|W| 2 + V tr m 2 + i 5 |*| 4 + 



(114) 
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in the rotating frame (for simplicity, we now use E instead of E'). In Sec. III.B, a physically motivated wave function 
served to evaluate Eq. (114), yielding Eq. ( [49| ) for the energy of a straight vortex displaced laterally from the trap 
axis. As noted previously, the assumption of a straight vortex restricted the analysis to a disk-shape condensate. 

To find the energy of a curved vortex, one can first find the condensate wave function ^ and then substitute it 
into the functional (114). For a curved vortex line, however, this approach is complicated. Instead, one can use 
Eq. (113) to find the vortex energy directly. As we know, the stationary Gross-Pitaevskii equation can be obtained by 
varying the energy functional (114). The dynamical equation ( |113[ ) is, in fact, the time dependent Gross-Pitaevskii 
equation, written in a way suitable to describe the vortex motion. Consequently, if we formally put V(r) = in 
Eq. (113) (namely omit the time derivatives), then the resulting stationary equation must be an extremum of the 
energy functional Ey associated with the presence of the vortex and considered as a functional of the vortex shape 
Ey = Ey (x(z), y(z)). An equivalent energy functional has the form (in the TF limit): 



r , , , 

Ey {x{z),y{z)) = — 



dz 



|*n|-\ I + + (//)- In 



2M g\q TF \ 4 n 



A ± V tI 



(115) 



where the prime denotes the derivative with respect to 



gives Eq. (113) with V(r) = 0, apart from terms of higher order xx' 2 ,xy' 2 



Variation of Eq. (115) with respect to x{z) and y{z) 
. Hence Eq. (115) provides an energy 



functional for the small deformations of a vortex about a straight configuration along the z axis (when the fourth-order 
terms in the displacement can be omitted) or for arbitrary displacements of a straight vortex. Note that Eq. ( |115| ) 
involves only a one-dimensional line integral instead of the three-dimensional expression in Eq. (114), which is a 
significant simplification. In scaled dimensionless units x — > R x x, etc., this energy functional becomes 



Ey (x(z),y(z)) = 27r^i?^ 2 n(0) / dz 



1 



y 2 - z 2 ) y/l + a(x>) 2 + (3(y') 2 In ( ^ ) - 



2/xft(l - x 2 - y 2 - z 2 f 



(116) 



where n(0) = fi/g is the density at the center of the vort ex-fr ee condensate, £ 2 = fi 2 /2M/i and the integration is 
restricted to the region 1 — x 2 — y 2 — z 2 > 0. Using Eq. (11£) one can obtain a simple expression for the angular 



momentum of the condensate in the presence of a curved vortex line: 



dEy 

on 
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(117) 



where N = 8TrR x R y R z n(0)/15 is the total number of particles in the condensate. 

The integration in Eq. ( |ll6j ) is particularly easy for a straight vortex and readily reproduces Eq. (p9[). An expansion 
for small lateral displacements yields Eqs. (|5l]) and (p3^) for il m and f2 c for a disk-shape TF condensate. In the more 
general case of arbitrary small displacements, Eq. ( |116[ ) can be expanded to second order in the amplitudes x and y 
and their derivatives. Use of the dynamical equations that lead to (100) and (|101|) gives the simple expression 



Ey (x(z),y(z)) = y^^"(0) 
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(118) 



The first term of Eq. (118) reproduces the value of Q c for a general TF condensate, and the second term becomes a 
sum over all normal modes of the form (B9f) 



Ey(x(z),y(z)) = ^-vR z en(0) 



5 h{u 2 x 



dz(l~z 2 )x n (z)y n (z), (119) 



where the orthogonality condition Eq. (104) eliminates the cross terms between different normal modes. If any of the 
normal modes is anomalous (namely with negative frequency) , then the system is unstable with respect to excitation 
of those modes. This analysis confirms the interpretation of f2 m as the applied rotation frequency at which the 
frequency of the last anomalous mode vanishes in the rotating frame. At this applied O the location of the vortex 
line along the z axis becomes a local minimum of energy. Note that this conclusion is wholly equivalent to that in 
Eq. ( |60| ) based on the Bogoliubov quasiparticles. 

One should note that for a cigar-shape condensate with R z > 2R± , there is an interval of angular velocity of trap 
rotation when Q c < Q < Q m . In this interval, the frequency of (at least) the lowest vortex mode remains negative, 
but penetration of a vortex into the condensate is energetically favorable. Under such a condition, the vortex line can 
lower its energy by undergoing a finite-amplitude deformation, and the ground state of the system corresponds to a 
curved vortex line displaced from the trap axis (see also [116]). 
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4- Precession and tilting of a straight vortex line in a nearly spherical TF condensate 



The preceding discussion of vortex dynamics in a three-dimensional confined condensate has focused on the small- 
amplitude displacements from equilibrium. In the special case of a spherical trap, however, the presence of a zero- 
frequency precessing mode (Sec. V.D.2) allows a more general analysis of the nonlinear dynamics, which is directly 
relevant to recent JILA experiments on the evolution of an initially straight vortex in a nearly spherical TF conden- 
sate |9^]. In practice, the trap deviates slightly from spherical with R x ^ R y ^ R z . 

For a spherical condensate, a motionless straight singly quantized vortex through the center of trap satisfies the 
general Eq. (|9^) for the velocity of a vortex line because the axis of the vortex i lies along VVtr- Let 

x = jxs, y = 7 y s, z = j z s (120) 

specify the axis of the vortex line, where s is the arc length measured from the trap center and (7^, 7j,, 7z) are 
the direction cosines relative to the principal axes of the anisotropic trap. For small anisotropy, the vortex remains 
approximately straight, but the direction cosines become time dependent. To first order in the anisotropy, the 



curvature k can be omitted in Eq. (93) and I^tfI can be approximated by the TF density for a spherical vortex-free 



^taffW-afl^,, (121) 



condensate with TF radius R. Standard perturbation theory yields the nonlinear dynamical equations 

lx = 7- 



This set of equations is familiar in classical mechanics as Euler's equations for the torque-free motion of a rigid 



body [131 -133 1, where they describe the motion of the angular- velocity vector as seen in the body-fixed frame. In the 



present context, this set of three coupled nonlinear equations has two first integrals 

= ^ (124) 

which verifies that the first-order anisotropy simply rotates the vortex axis and 

ulll + u 2 y ll + ^lll = const, (125) 

which is the condition of energy conservation. 

The simplest situation is an axisymmetric trap with lj x — lj v = u± , in which case the v ortex line precesses uniformly 



about the z axis (the symmetry axis) at a fixed polar angle arccos7 2: (0) at a frequency [117| 
5h(w* - u)\) (1MR\ 5h ( 1 1 \ (IMR 



4, *<°> ) - 2M U " Rl ) ^ ln {— ' ' (126) 



where the numerical factor 1.96 inside the logarithm is the same as that discussed below Eq. (107). For positive (neg- 
ative) lo, the precession is counter-clockwise (clockwise). Recent experiments at JILA have observed two recurrences 
of such precessional motion in a slightly flattened trap with oj z — lu±_ w 0.1uj z and a polar tipping angle of 45° from 



the z axis. In this case, Eq. (126) predicts u/2ir w 0.33±0.03 Hz, in an agreement with the observed value 0.25±0.02 
Hz @. 

More generally, for an anisotropic trap (with uo x > uj y > u) z ), the vortex executes closed trajectories (see Fig. |l3|). For 
initial positions close to the x and z axes (the smallest and largest TF radii), the motion is "stable," remaining nearby, 
but small-amplitude motion about an initial position close to the y axis (the intermediate TF radius) yields imaginary 
frequencies. Thus such trajectories deviate far from the initial neighborhood, even though they eventually return 



(this periodic behavior is familiar from the corresponding solutions of the Euler equations [131-133]). Reference [117 



gives explicit solutions for the resulting dynamical motion of a nearly straight vortex in a totally anisotropic trap. 



34 



FIG. 13. Typical trajectories of the end of a straight vortex line (that passes through the condensate center) during its 
motion in a slightly nonspherical trap with R x < R y < R z . 



VI. EFFECT OF THERMAL QUASIPARTICLES, VORTEX LIFETIME AND DISSIPATION 

In previous sections we considered a Bose condensate within the Bogoliubov approximation, which omits the effect 
of thermal quasiparticles. At finite temperatures, however, these noncondensate atoms can modify the frequencies of 
the vortex modes and dissipate energy. 



A. Bogoliubov and Hartree-Fock-Bogoliubov Theories 



Let us consider a condensate in thermal equilibrium at temperature T. Within the Hartree-Fock Bogoliubov (HFB) 
theory, the condensate wave function s atisfies the following generalized Gross-Pitaevskii equation (in a frame 
rotating with the angular velocity ttz) [134]: 



~ 2M + VtI + 9m + 2gpir) ~ ^ + iMld * I * + aA ^* = °' 



(127) 



where <fi is the azimuthal angle in cylindrical polar coordinates, p(r) is the density of the noncondensed gas and A(r) 
is the anomalous average of two Bose field operators describing the noncondensate (as in Sec. IV.A, ip = \P + <fi is the 
quantum field operator, with A = (0</>) and p — ((f)* </>)). The collective excitation energies E of the system are the 
eigenvalues of the generalized Bogoliubov equations for the coupled amplitudes u(r) and v(r) 



2M 



X7 2 + V tI + 2g\y\ 2 + 2gp(r)~ f i(n) 



-g [A*(r) + ** 2 ] 



-g [A(r) + * 2 ] 
—iMlds 



= E 



u 



(128) 



Eq. ( 128 ) is valid at least for temperatures much less than the chemical potential p w hen resonant contributions 
(the so-called Szepfalusy-Kondor processes) to the self-energies are not substantial [135]. In addition, we have self- 
consistency relations for the noncondensate density p(r) 



P(r) = E 



exp(E n /k B T) - 1 



K(r)| 5 



(129) 
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and for the anomalous average A(r) 



A(r) = -]T 



2u n r )v* r 

+ u n {r)v n {r) 

exp(E n k B T) - 1 



(130) 



where n denotes quantum numbers specifying the excited states with energies E n (n — 0, 1, 2, • • • ). The eigenfunctions 
u n (r) and v m (r) satisfy the normalization condition: 



[u* (r)u m (r) - u*(r)v m (r)] rfr = 5 n 



(131) 



Equations ( 127 )-( 131 ) constitute a complete set of the self-consistent equations for the HFB theory. Within this 
theory, the quasiparticle eigenvalues E n in Eqs. (128)- (130) must be positive because the condensate is defined to have 
zero energy. Thus a negative eigenvalue means a failure of the self-consistency and the associated thermal equilibrium 
of the system. If p(r) and A(r) are set to zero, we recover the Bogoliubov theory. If we set only A(r) = 0, we obtain 
the Popov approximation. For a vortex-free condensate in the low-temperature limit, the Popov and Bogoliubov 



theories give identical excitation spectra |138|. The excitation spectrum in the HFB theory has an unphysical gap 
because it does not treat all condensate-condensate interactions consistently [ 134 1 . Gapless modifications of the HFB 
theory, the so-called Gl and G2 approximations, are discussed in [ 139, 140 1. Normally, the zero-temperature limit of 
the Popov, Gl and G2 theories should be the Bogoliubov theory (which does not take into account noncondensate 
atoms). For a non rotating condensate with a vortex, however, this is not the case because vortex is unstable. 

Within the Bogoliubov theory, an isolated vortex in a nonrotating harmonic trap has at least one normal mode 
with negative energy. Let us apply the HFB theory for a condensate with a vortex. To find a self-consistent solution 
for t he lo west eigenvalue at low temperatures, one can use a perturbation method analogous to those developed in 
Ref. [ 114 1 . We consider a condensate in an axisymmetric trap that rotates with an angular velocity £1 around the z 
axis. We assume that the condensate contains a singly quantized vortex along the z axis. For simplicity we consider 
a disk-shape condensate, so one can omit vortex curvature in investigating the lowest normal mode. The condensate 
wave function has the form ^ = e^l^l, with A = e 2l( ^|A|, and we can rewrite the generalized Bogoliubov equations 
as: 



where 



J^v* + 1 -Muy + 2 9 \* f- tt{ n))(l I] ry _. ma 



(132) 



(133) 



xy plane with the same chemical potential; its excitations obey the equation 



and V includes the remaining part of Eq. (125). Here, is the wave function for an unbounded condensate in the 

(134) 



H 



u 
-vo 



Equation (134) has an exact pair of solutions (see Ref. [114]) with positive norm and energy 

E = Ml. (135) 
Let us now make the following assumption: E -C kgT -C Ei, E2, ■ ■ •, where E is the energy of the lowest normal 



mode, which can depend on T. Then the term with n = gives the main contribution in the sum in Eqs. (|129j), (130) 
and we obtain: 



p(r)«^[K(r)| 2 + Mr)| 2 



(136) 



A(r)«-^u (rK(r) (137) 



For a singly quantized vortex one can derive the expression 
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where I 2 f=a 16v / 2~7'"A t3 ^ 2 /3<7 t '-'zV^ : is a normalization integral, and 

|A(r x =0,*)| «0. 

In first-order perturbation theory, the lowest energy eigenvalue Eq is defined by the equation: 

e = m + E a - 1 



(138) 



(139) 



(140) 



where E a = (3H 2 w 2 _/4:fi) In (R±/Q and 7 = 0.077i?i/A^ 4 are positive with N = ^7rn(0)i? z i?^ the total number of 
particles in the condensate, and fi can be taken as the chemical potential for a nonrotating trap. Eq. (14C) has two 
solutions, one with positive energy and one with negative energy that reproduces the previous anomalous mode with 
Eq = Ml — E a as T — > 0. The negative solution can be formally omitted, satisfying the requirement of self-consistency. 
The positive solution has the form: 



1 



En = 



For nonrotating trap (O = 0), we find 
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If T — ► 0, we obtain Eq ps jkgT, so that Eq is proportional to T in the low-temperature limit. In fact our method 
generalizes the Beliaev theory ]136| for the vortex state. Recently Pitaevskii and Stringari actually generalized the 
Beliaev approach (in the density-phase representation) for the trapped Thomas-Fermi condensate [ 137 1 . 

Virtancn, Simula and Salomaa made numerical calculations of vortex normal modes at finite T within the Popov, 
Gl and G2 approximations and demonstrated that for a singly quantized vortex there is a self-consistent solution 
with only positive frequencies in the limit T — > [ 141 1 . Their lowest energy solution corresponds to our Eq. (142). 
The vortex mode ( |142[ ) arises from the presence of quasiparticles (an external pinning potential can also result in such 
motion [ 142 1 ) . At low temperatures, the quasiparticles are mostly localized in the vortex-core region and provide an 
extra repulsive potential [the term 2gp(v) in Eq. (128)] that affects the elementary excitations. At T = 0, the residual 
localized noncondensate fraction arises from the interaction between particles; this result follows from Eq. ( |138| ) if we 
take Eq oc T at low temperatures. The additional potential has a peak at the vortex core and the vortex line precesses 
around the quasiparticle potential center with a positive excitation energy. 

However, this does not mean that quasiparticles stabilize the vortex in a trap. The physics of the problem is 
the following. At any moment during the vortex motion, quasiparticles fill the vortex core (the relaxation time of 
quasiparticles is much less than the period of the vortex precession). The vortex line participates in two motions: 
first, the vortex precesses around the trap center with the frequency fku a = —E a < (O = 0). The trap potential is 
responsible for this unstable mode. The quasiparticles are localized in the vortex core and move together with the 
vortex; their presence simply slightly changes the chemical potential and slightly decreases the normal mode frequency. 
In second motion, the vortex line moves around the center of mass of the quasiparticles in a locally uniform condensate 
(in xy plane). The amplitude of this motion is less than £ and the frequency can be found from Eq. (142) in the limit 



— > 00 or E a 



0: 



(143) 



where Uq is the density of the vortex free condensate at the vortex location (in the plane z = 0). For JILA parameters 
7 w 0.3, E a w 1.58Hz, then for T = 0.8T C we obtain lux ~ 13-6 Hz. If this mode is thermally excited, its amplitude is 
given by 



A = e 



Go 



1/2 



fk B T\ 
\ Hujt J 



1/2 



(144) 
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where a is the scattering length. For parameters of JILA experiments A = 0.16£. Taking into account lot <x VT we 
obtain the following temperature dependence A oc T 1 / 4 . It is interesting to note that the thermal mode ( |143[ ) exists 
only in 3D condensate; in the limit R z — oo both the mode frequency and the amplitude go to zero. 

Recent measurements of the lowest vortex modes in the JILA experiments are in a good quantitative agreement 



with solutions of the time-dependent Gross-Pitaevskii equation |3^,|83,117]. The JILA experiments measure, in fact, 
not only the absolute value, but also the sign of the lowest vortex mode. The negative value of the anomalous-mode 
frequency means that the vortex precesses in the same direction as the superfluid flow around the vortex core, which 



is seen in the experiments. An experimental observation of the thermal mode (143) could be next challenging problem 
of future investigations. 

B. Dissipation and Vortex Lifetimes 

It is valuable to consider dissipation and its role in the vortex lifetime. In a nonrotating trap, the ground state of 
the system is a vortex-free condensate, so that a condensate with a vortex necessarily constitutes an excited state. In 
the absence of dissipation, however, the vortex line moves along trajectories of constant energy, remaining inside the 
condensate. The condensate with a vortex will be unstable only if there is a mechanism to transfer the system to the 
lower-energy vortex-free state [BTj . The dissipative dynamics of a straight vortex due to its interaction with the thermal 



cloud in a trapped Bose-condensed gas was discussed by Fedichev and Shlyapnikov [ 143 1 . If the vortex line moves with 
respect to the normal component, scattering of elementary excitations by the vortex produces a friction force, like 
that in superfluid 4 He (see Ch. 3 of Ref. |l4|]). Such a mechanism can transfer energy and momentum to the thermal 
cloud. The friction force F can be decomposed into longitudinal and transverse components: F = Du — D'(u x h), 
where u is the velocity of the vortex line with respect to the normal component, D and D 1 are the longitudinal 
and transverse friction coefficients, respectively, and h is a local tangent vector to the vortex line. The transverse 
friction coefficient is independent of the scattering amplitude and is given by the universal expression D' = hp n /M , 



where p n is the local mass density of the normal component 144]. The longitudinal friction coefficient depends on 
the scattering process. In the limit fc^T ^> /i, one can treat the elementary excitations as single particles, with 
the result that p n ~ 0.1M 5 / 2 T 3 / 2 //i 3 and the longitudinal friction coefficient is proportional to the temperature: 
D ps (na 3 ) 1 / 2 ! 1 /^, where n = |^| 2 is the superfluid density for the vortex-free condensate and a is the s-wave 
scattering length |143[ |. 

In the presence of dissipation, the vortex line moves toward a (local) minimum of the energy. In a nonrotating 
condensate, an off-center vortex precesses around the trap center and is expected to spiral out to the condensate 
boundary due to the dissipation. Once the vortex reaches the boundary, it presumably decays by emitting phonons 
and single-particle excitations. The radial motion of the vortex is governed by the longitudinal friction coefficient: 
v r Du/hn ~ (na 3 ) 1 / 2 Tw//i <C u, wher e u is the precessional speed. Using this expression, one can estimate the 



characteristic lifetime of the vortex state [143]. At present, no dissipation of the moving vortex has been observed in 
the JILA experiments [fS7| . A characteristic decay time for the dissipative mechanism of Fedichev and Shlyapnikov in 
the JILA conditions is significantly larger than the life-time of the condensate. The temperature and density are too 
small to see the dissipation. 

Another factor that can influence the vortex lifetime is the possibility that a moving vortex can emit phonons. It 



is known that a moving vortex in an infinite compressible fluid emits phonons, leading to a slow loss of energy [145 



Recently, Lundh and Ao [123] studied the radiation of sound from a moving vortex in an infinite, uniform system. 
A homogeneous two-dimensional superfluid described by a nonlinear Schrodinger equation is equivalent to (2+1)- 
dimensional electrodynamics, with vortices playing the role of charges and sound corresponding to electromagnetic 



radiation [146,147]. Thus, a vortex moving on a circular trajectory in an infinite superfluid radiates sound waves, 
which are analogous to the cyclotron radiation of an electrical charge moving along a circular orbit. The power 
radiated by a vortex with unit length executing circular motion with frequency to at a radius ro is given by the 



following Poynting vector [123] 



P-I^S, (145) 

where Q = is the "vortex charge," n is the uniform superfluid density, and c s = y/ /i/M is the velocity 

of sound. 

In a nonuniform system, such as a two-dimensional or a disk-shape axisymmetric trapped condensate, an off-center 
vortex performs a circular motion around the symmetry axis. If such motion excites sound waves (radiates energy), 
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the vortex will move outward toward regions of lower potential energy, until it eventually escapes from the cloud. In 
a trapped condensate, however, the excitations all remain confined within the condensate, and no phonon radiation 
is expected. In particular, the wavelength A of sound that would be emitted exceeds the size R of the condensate. 
Indeed, A ~ 2ttc s /uo and the precession frequency of the straight vortex is of the order of u> ~ hh\(R/£)/MR 2 \ as a 
result, X/R ~ (R/Q ha(-R/0 ^ L an d the "cyclotron" radiation is prohibited. 

Finally, let us discuss how vortex generation affects the dissipation in superfluids. One classic manifestation of 
superfluidity is that objects traveling below a critical velocity propagate through a superfluid without dissipation. 



According to the Landau criterion 1 105 ] , which relics on the use of Galilean invariance, the critical velocity is vl = 
Tain[E(p)/p\, where E(j>) is the energy of an elementary excitation with momentum p. For a homogeneous Bosc 
condensate, the Bogoliubov spectrum implies a Landau critical velocity equal to the speed of sound vl — c s . The 
Landau critical velocity can usually be observed only by moving microscopic particles through the superfluid. Such 



motion of microscopic impurities through a trapped gaseous Bose condensate was studied recently in [ 148 1 . As the 
impurities traverse the condensate, they dissipate energy by colliding with the stationary condensate and radiating 
phonons. When the impurity velocity was reduced below the speed of sound, however, the collision probability 
decreased dramatically, providing evidence for superfluidity in the condensate. 

If a macroscopic object moves through the condensate, dissipation can occur due to turbulence and vortex formation 
in the superfluid, even if the object's velocity is much lower than the Landau critical velocity. Recently, dissipation in 



a Bose-Einstein condensed gas was studied by moving a blue-detuned laser beam through the condensate |149 , |l50 |. 
The laser beam repels atoms from its focus and creates a moving macroscopic "hole" in the condensate. The observed 
heating of the system agrees with the prediction of dissipation when the flow field becomes locally supersonic. Nu- 
merical simulations of the nonlinear Schrodinger equation were used to study the flow field around an object moving 



through a homogeneous condensate [28 151-153.122]. When the object moves faster than a critical velocity v c , these 
studies show that the superfluid flow becomes unstable against the formation of quantized vortex lines, which gives 
rise to a new dissipative regime. Pairs of vortices with opposite circulation are generated at opposite sides of the 
object. The rate of the energy transfer to the condensate by the moving object increases significantly above this 
critical velocity for vortex formation. The heating rate can be expressed as dE/dt — E pa x T f s , where -E pa i r is the 
energy of a vortex pair and f s is the shedding frequency. The rate of vortex-pair shedding f s is proportional to v — v c 
and thus larger when the speed of sound is lower. 

Other simulations of the GP equation have demonstrated that vortex-antivortex pairs or vortex half-rings can be 



generated by superflow around a stationary obstacle [152,154,155 2q] or through a small aperture |156|. One might 



expect similar excitations in a rotating condensate. In addition, vortex half-rings can be nucleated at the condensate 
surface when the local tangential velocity exceeds a critical value. 

VII. VORTEX STATES IN MIXTURES AND SPINOR CONDENSATES 



The advent of multicomponcnt BECs [ 157 - 159 has provided many new possibilities for quantum- mechanical state 
engineering. Since there is no intrinsic difficulty in loading and cooling more than one alkali element in the same 
trap, interpenetrating superfluids can now be realized experimentally. Binary mixtures of condensates can consist of 
different alkalis, or different isotopes, or different hyperfine states of the same alkali atom. Such binary mixtures of 
Bose condensates have a great variety of ground states and vortex structures that are experimentally accessible by 
varying the relative particle numbers of different alkalis In particular, one can move continuously from regimes 
of interpenetrating superfluids to those with separated phases. Many alkali binary mixtures contain a coexistence 



region, which is the analog of He- He interpenetrating superfluids in ultralow-temperature physics [16C 



A. Basic Phenomena 

Most experiments on Bose-Einstcin condensation of atomic gases of 87 Rb §, 7 Li ||, and 23 Na § have used magnetic 
traps to condense atoms with a hyperfine spin F = 2 (or F = 1). Such a condensate of spin-F bosons constitutes a 
spinor field 

$ ro (r,t))=C m (r,t)¥(r,t), (146) 

where ip m is the field operator, m labels F z (where — F < m < F), \& is a scalar, and C, m is a normalized spinor. 
In magnetic traps, the spins of the alkali atoms are frozen and maximally aligned with the local magnetic field B 
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||. As a result, £ is given by the eigenvalue equation B • F£ = F£, where F is the hyperfine spin operator and B 
is a unit vector along B. The dynamics of (ip m ) is therefore completely specified by the scalar field 'J, as in He. 
Thus, even though the alkali atoms carry a spin, they behave in magnetic traps like scalar particles. In contrast to 



the scalar field, however, the spinor field in Eq. ( J146[ ) possesses a local spin-gauge symmetry: a local gauge change 
exp[ix(r, t)] of (V> m ) can be undone by a local spin rotation exp[— i(x/F)fi(r, t) ■ F]. Because of this symmetry, the 
effective Hamiltonian of the scalar field \? is not that of 4 He, but that of a neutral superfluid in a velocity field u s . 
The velocity (or gauge field) u s is a direct reflection of the spin-gauge symmetry and it is given by 

u. = -^CM. (147) 



The velocity u s can be calculated from the vorticity f2 s of u s , which satisfies the Mermin-Ho relation 161 162], 



Equation (148) shows that the spatial variations of B necessary to produce the trapping potential will inevitably 
generate a nonvanishing superfluid velocity u s — (2h/M) (1 — B z /B) V [a.Ycta.n{B y / B x )\ 0. If Bo = Bqz is the 
magnetic field at the center of an axisymmetric harmonic trap and ujq is the maximum trap frequency, then the 
spin-gauge effect generates the following constant effective "rotation" fl s around the z axis || : 

f2 = . hlOn 



(149) 



where fiB is the Bohr magneton. The superfluid velocity u s splits the degeneracy of the harmonic energy levels, 
breaks the inversion symmetry of the vortex- nucleation angular velocity Q c , and can produce vortex ground states in 
the absence of external rotation if f2 s > il c Q. In current experiments, the spin-gauge effect is small; for example, if 
uio = 10 Hz and Bq = 1 G, we obtain £I s /ujo ~ 10~ 5 . In oblate traps with lo z ^> uj±, however, the spin-gauge effect 
can be significant (f2 s could be comparable with to± for large enough values of ui z ). 



Recently, the MIT group has succeeded in trapping a 23 Na Bose condensate by purely optical means |l58| , |l59 |. 
In contrast to a magnetic trap, the spins of the alkali atoms in such an optical trap arc essentially free, so that the 
spinor nature of the alkali Bose condensate can be fully realized. Specifically, 23 Na atoms possess a hyperfine spin, 
with F — 1 in the lower multiplet. All three possible projections of the hyperfine spin can be optically trapped 
simultaneously. Thus the condensate is described by a spin-1 spinor. The internal vortex structure of a trapped 



spin-1 BEC was investigated in Ref. [ 163 1 . Such vortices and their stability were also discussed in [20 164 1. In an 
optical trap, the ground state of spin-1 bosons such as 23 Na, 39 K, and 87 Rb can be cither ferromagnetic or "polar," 
depending on the scattering lengths in different angular momentum channels pp[ |. The ferromagnetic state also has 
coreless (or skyrmion) vortices, like textures found in superfluid 3 He-A. Because of the wide range of hyperfine spins 
of different alkalis, the optical trap has provided great opportunities to study different spin textures in dilute quantum 
gases of atoms with large spins. This is a fruitful subject for future experiments. 

Although most of the theoretical effort has concentrated on single-condensate systems, the first experimental real- 



ization of BEC vortices was achieved with a two-species 87 Rb condensate |33|], following the proposal of Ref. [165|. 
Several other proposals have been made for the dynamical production of a vortex using the internal structure of atoms 
[ 166 169 1 . The spin-exchange scattering rate is suppressed for 87 Rb, which makes possible the study of magnetically 



trapped multicomponent condensates of these atoms. The two species correspond to two different hyperfine energy 
levels of 87 Rb, denoted |1) and |2); they are separated by the ground-state hyperfine splitting. Since the scattering 
lengths are different, both states are not equivalent. Typically, the |1) = \F = l,m = —1) state is trapped and cooled 
to the condensation point. Once the atoms in |1) have formed the condensate ground state, a two-photon microwave 
field is applied, inducing transitions between the |1) state and the |2) = \F = 2, to = 1) state [B3| . As a result, the 
atoms cycle coherently between the two hyperfine levels with an effective Rabi frequency f2 c ff Jl70[ . Two parameters 
characterize the coupling: the detuning and the power. The detuning S denotes the mismatch of the frequency of 
the coupling electromagnetic field to the frequency difference between the two internal atomic states. The power is 
characterized by the Rabi frequency f2; it is the rate at which population would oscillate between the two states if 5 
were zero. When 5 is larger than f2, the population oscillations occur at the effective Rabi frequency £l c g — yjfl 2 + S 2 , 
which obviously exceeds Q. 

In principle, both states could be cooled simultaneously, so that the condensate forms in a mixture of states. In 
practice, however, the typical lifetime of atoms in the |2) state is about 1 s due to inelastic spin-exchange collisions, 
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which makes it very difficult to achieve runaway evaporation for this state. In contrast, atoms in the |1) state have 
a much longer lifetime of about 75 s Q. The advantage of using the \F — l,m = —1) and |F = 2, to = 1) states 
is that their magnetic moments are nearly the same, so that they can be simultaneously confined in identical and 
fully overlapping magnetic trap potentials. Unlike the more familiar single-component superfluids [see a discussion 
after Eq. J153Q ] , where the topological constraints make it difficult to implant a vortex within an existing condensate 
in a controlled manner, the coupled two-component condensate has a different order parameter and hence different 
topological constraints. Indeed, the coupled two-component system allows the direct creation of a |2) (or |1)) state 



wave function having a wide variety of shapes out of a |1) (or |2)) ground-state wave function [ 165 1 . 

For example, to form a vortex in the two-component system, one should impose a perturbation H\ that couples 
the ground state of the system to the vortex state (namely, the matrix element of the perturbation operator between 
these two states must be nonzero). The time-dependent GP equation describing the driven, two-component condensate 



is [165 



ih!L f^A - (H Q + U 11 \y 1 \ 2 + U 12 \y 2 \ 2 + H 1 + hS/2 m/2 nnm 

dt\*2j \ fin/2 H + U 21 \^ 1 \ 2 + U 22 \^ 2 \ 2 -H 1 -M/2J\^ 2 J' U ° Uj 

where Ho = — (fi 2 V 2 /2M) + \MlJq{t\ + z 2 ) for a spherical trap, M is the atomic mass, u>q is the trap frequency, 
Uij = 47rfi. 2 ajj/M, with Oij the s-wave scattering lengths for binary collisions between constituents i and j. Williams 
and Holland considered the perturbation Hi in the following form |165H : 



Hi = «[/(r) cos(wi) + g(r) sm{wt% (151) 

where k is a coupling coefficient and /(r) and g(r) are prefactors that depend on r. The explicit form of Hi determines 
the symmetry of the quantum state being prepared, so that general / and g can serve to prepare a macroscopic 
quantum state of arbitrary symmetry. To create a vortex state with one unit of angular momentum, one can take 
k — MujqPq, /(r) = x and g(r) = y in Cartesian coordinates. This form of perturbation effectively confines the two 
hyperfine states in separate axially symmetric harmonic-oscillator potentials with the same trap frequency u)q. The 
trap centers are spatially offset in the xy plane by a distance po (from the center) and rotate about the symmetry axis 
at an angular velocity ui. To achieve this configuration experimentally, Ref. [^3| shone a laser beam into the trap along 
the z axis so that the cloud sits in the middle of the Gaussian beam waist where the gradient of the beam intensity 
is approximately linear (see Fig. |l4|a). This arrangement produces a constant force on the atoms. If the frequency 
of the laser beam is tuned between the two hyperfine states, the optical dipole force acts in opposite directions for 
each state, displacing the trap centers for each state. When the beam rotates around the condensate at the angular 
velocity u), we obtain the desired result. 

To create a vortex, the angular velocity u> should be close to the value at which a resonant transfer of population 
from the nonrotating condensate into the vortex state takes place. Consider the frame co-rotating with the trap 
centers at an angular frequency w. In this frame, the energy of the vortex with one unit of angular momentum is 
shifted by hui relative to its value in the laboratory frame. When this energy shift compensates for both the energy 
mismatch HS of the internal coupling field and the small chemical potential difference between the vortex and the 
nonrotating condensate, resonant transfer of population takes place (see Fig. |l4|b). It is obvious that if we change the 
sign of detuning S while keeping the trap rotation fixed, a vortex will be created with opposite circulation. Vortices 
with opposite circulations experience opposite energy shifts in transforming to the rotating frame and therefore require 
opposite signs of detuning in order to achieve the resonant coupling. 

In practice, ui S> loq and 5 Q. The first inequality allows the vortex to be generated rapidly. The main problem 
with a slow drive (when lu m lu ) is that the time scale for coupling to the vortex state is very long, on the order of 
seconds in a trap with luq — 10 Hz. The weak-coupling limit, given by the second inequality, allows the resonance 
condition uj w 5 to select energetically the desired state with high fidelity. 



Figure |l5j shows the results of a numerical integration of Eq. (150) in two dimensions [u> z = 0), with the condensate 
initially in the nonrotating ground state and in the internal state |1) [165]. The coupling drive is turned on at time 
t = 0, and is turned off at time t = t s by setting both fl and po to zero. The top and the bottom graphs show the 
fractional population and the angular momentum per atom of the |2) state as a function of time. 

The small-amplitude rapid oscillations on the top graph correspond to the cycling between internal levels at the 
effective Rabi frequency fi e ff • The gradual rise of this line reflects coupling from the ground state to the vortex mode 



caused by the drive Hi in Eq. (150). Once during each Rabi cycle, the angular momentum approaches unity (bottom 
graph), and, at that time, the |2) state wave function approaches a pure vortex mode. By turning off the coupling 
at a precise time t = t s on a given Rabi cycle, the \2) state can be prepared to have unit angular momentum. The 
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FIG. 14. (a) A basic schematic illustration of the technique used to create a vortex. An off-resonant laser provides a rotating 
force on the atoms across the condensate as a microwave drive of detuning 8 is applied, (b) A level diagram showing the 
microwave transition to very near the \2) state, and the modulation due to the laser rotation frequency that couples only to 
the angular momentum I = 1 state when uj ~ 8. 
(Taken from Ref. @) 
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FIG. 15. Dynamical evolution that can create a vortex. The top graph shows the fractional population of atoms in the |2) 
internal state. The bottom graph shows the angular momentum of the 2) state, in units of Planck's constant h. The inset 
shows the amplitude of population transfer to the vortex as a function of the trap rotation frequency uj, with A = f2 e fi — 
The various parameters used in the calculation are: uja = 10 Hz, 8 = 200 Hz, uj — 205.4 Hz, N — 8 x 10 5 atoms, M is the mass 
of the 87 Rb atom, for simulations the values of scattering lengths are taken to be an = 022 = ai2 = 5.5 nm, and for t < t s , 
Q = 50 Hz and po = 1.7/im. Reprinted by permission from Nature 401, 568, (1999), ©1999 Macmillan Magazines Ltd. 
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maximum possible population transfer to the vortex state using this scheme obeys a Lorentzian response curve as u) 
is varied near f2 e ff, exhibiting a narrow resonance. This situation is shown in the inset of Fig. [l5], where A = fl c g — to. 

In an experiment, it is possible put the initial condensate into either the |1) or |2) state, and then make a vortex 
in the |2) or |1) state, respectively. The evolution of the vortex can be watched over time scales from milliseconds to 
seconds. In Ref. [ p3| , the vortex was found to be stable in only one of the two possible configurations corresponding 
to the vortex in the |1) state, which is the one with the larger scattering length (with the |2) state in the core). The 
other possibility (the vortex in the |2) state, which is the one with the lowest self-interaction coefficient) produces an 
instability. 



B. Stability Theory 



We use the following notation for the states: (1, 0) for the state with the vortex in |1) and (0, 1) for the state with 
the vortex in |2). In the JILA experiment fjS], the number of particles is the same for each component (Ni = iV 2 = N) 
but, in general, one could consider any ratio between the populations of the different levels. The scattering lengths for 
binary collisions depend on the internal hyperfine level of the atom. For 87 Rb the values of scattering lengths are nearly 
degenerate and in the proportion an : ai 2 : a 22 = 1-00 : 0.97 : 0.94 [171]. Because of the relation Un > U\i > U21, 
the experiment is performed in a regime in which the first component separates from the second one. Consequently, a 
favored configuration has the first component spread over the largest part of the space. Numerical simulations show 
that in the equal population case, Ni = A 2 = N, and for arbitrary nonlinearities, the stationary states (1, 0) is stable 
while the other state (0, 1) is unstable. 



The origin of the instability of the state (0, 1) is purely dynamical [172] and can be understood within the framework 
of mean-field theories for the double-condensate system without dissipation. Actually, the instability mechanism does 
not lead to expulsion of the vortex from the condensate, but to periodic transfer of the phase singularity from one 
species to the other. To study the vortex stability, one can start from a pair of coupled Gross Pitaevskii equations for 
the condensate wave functions of each species 



ih— = 
ot 



ft 2 V 2 
' 2M 



^i + C/ii|*i| 2 + C/i 2 |* 2 | 



(152) 
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where V\ and V2 are trap potentials for the condensate components. These equations are a particular case of Eq. 
( |150| ) when the drive is turned off (Hi, Q, 6 = 0). Equations ( |152| ) and (153) conserve the number of particles in 
each hyperfine level. However, the angular momentum of each component is no longer a conserved quantity, and 
the topological charge of each species can change through the time evolution. Instead, what is conserved is the total 
angular momentum of the system 



ih J (Pr^ld^i+ih J d 3 r*2^*2 



(154) 



As in the JILA experiments, we assume that both potentials are spherically symmetric and have the form Vi(r) = 
V 2 (r) = \Muj^(r\ + z 2 ). For stationary configurations in which each component has a well-defined value of the 



angular momentum, the time and angular dependence are factored out 



*i(r_L,Z, 



= e - i ^ h e i ' !i Ui(r±,z), 



(155) 



with i — 1,2. We focus on three particular configurations, which are the lowest energy states with vorticity (qi, g 2 ) = 
(0, 0), (1, 0), (0, 1). They correspond to the ground state of the double condensate, and to the single vortex states for 
the |1) and |2) species, respectively. 

Linear stability analysis of the three states gives the following results [172|. For the (0,0) state, the frequencies of 
all normal modes are positive, as expected for the ground state of the system. Among the normal modes of the (1, 0) 
family, there is a negative eigenvalue, which means that there is a path in the configuration space along which the 
energy decreases (this is just the analog of the anomalous mode in the one-component system with a vortex). This 
path belongs to a perturbation that takes the vortex out of the condensate. As in the case of a single-component 
condensate, however, the lifetime of the vortex state is only limited by the presence of dissipation (without dissipation, 
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FIG. 16. 

length d — y^h/Muio 



Evolution of the position of the phase singularity in the xy plane. Lengths are given in units of the trap characteristic 
(a) Phase singularity in |1), (b) phase singularity in |2). 



(Taken from Ref. [173|). 



the configuration is dynamically stable). Finally, in the (0, 1) family there are normal modes with complex frequencies. 
The shape of the unstable modes is similar to the energy-decreasing modes of the (1,0) family — that is, they are 
perturbations that push the vortex out of both clouds. The imaginary part of the eigenvalues implies that vortices 
with unit charge in |2) are unstable under a generic perturbation of the initial data, whereas those in |1) can be 
long-lived. This conclusion is consistent with the JILA experiments, where a vortex in the |2) species was found to 
be unstable 

Numerical simulations of the vortex behavior for large perturbations show that the linearly stable state (1,0) is 
robust and survives under a wide range of perturbations, suffering at most a precession of the vortex core plus 
changes of the shapes of both components [172]. This behavior arises in both two- and three-dimensional simulations. 
In contrast, the unstable configuration (0, 1) develops a recurrent dynamics. In the first stage, the first component 
and the vortex oscillate synchronously (the hole in |2) pins the peak of |1)). These oscillations grow in amplitude, and 
the vortex spirals out. Finally the first component develops a tail and later a hole which traps the second component. 
The hole is a vortex that has been transferred from |2) to |1). Though not completely periodic, this mechanism 
exhibits some recurrence, and the vortex eventually returns to |2). The preceding behavior persists even for strong 
perturbations in a two-dimensional condensate. However, for large perturbations of a three-dimensional condensate, 
the dynamics may lead to a turbulent behavior [172|. 

In Fig. [l6|, it is shown how a small initial perturbation makes the phase singularity in |2) spiral out of the system 
while a phase singularity appears in |1) and occupies the center of the atomic cloud. This dynamics is recurrent. 

The preceding results are valid for the equal population case, iVi = A^. For any ratio of the populations N1/N2 
and any values of the nonlinear coefficients Uij, the stability conditions are the following |173|: The configuration 
(1,0) is stable if 




> 1 



an 
ai2 



(156) 



For 87 Rb, the inequality ( p_56[) i s always satisfied, which proves that the configuration with a vortex in |1) is always 
linearly stable, as found in p3|. Note that the stability properties do not depend on the total number of particles 
but only on the ratio between the populations. 

The stability condition of the configuration (0, 1) is 
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This inequality fails for a certain range of Ni/N%. For the case of 87 Rb the unstable range is N1/N2 £ [0.73, 1.49], 
which means that certain choices of the population imbalance allow stabilization of the vortex in |2). These results 
predict the possibility of stable vortex states for various multiple-condensate systems [173]. 

Energetic considerations show that the extra degree of freedom associated with the second component allows a more 
intricate structure for the free-energy surface. As a result, in a two-component system, it is possible to achieve a local 
minimum in the free energy at the center of the trap 124]. The presence of such a minimum implies the existence of 
a region of energetic stability where the vortex cannot escape and might generate a persistent current. 



VIII. CONCLUSIONS AND OUTLOOK 



In this paper, we have provided an introductory description of vortices in trapped Bose condensed gases. The 
main conclusion of our analysis is that the vortex dynamics in such systems is well described by the time-dependent 
Gross-Pitaevskii equation (at least for low temperatures). The nonuniform nature of the condensate results in the 
appearance of anomalous vortex mode(s) with negative frequency and positive norm. Trap rotation shifts the normal- 
mode frequencies and can stabilize the vortex state. To date, experimental measurements of vortex dynamics and 
other properties of vortex states are in a good quantitative agreement with theoretical predictions based on solutions 
of the GP equation. Deviations from the mean-field predictions could arise when the gas parameter n\a\ 3 is not very 
small (semiclassical corrections to the mean-field approximation were calculated in [174 1 ) or from "mesoscopic" effects 
associated with the finite systems. However, there is no experimental evidence for these effects so far. 

We have been able to cover only part of the existing literature on vortices in trapped condensates. Among important 
issues that we have not discussed are: different methods of vortices generation and detection, kinetics of vortex 
nucleation and decay, vortices in BECs with attractive interactions and in Fermi condensates, other defects in BECs 
(solitons, instantons, vortex solitons, skyrmions, wave-function and spin monopoles). 

In the case of superfluid helium, vortex nucleation is associated with pinning of vortex lines at the walls of the 
container. Trapped condensates have no rough surfaces, and the nucleation process of quantized vorticity has a 
different origin [175, 116, 176 1. An important question in vortex nucleation is the role of the thermal component and 
transverse anisotropy of magnetic traps [177,178]. 

The literature of the past few years contains many different proposals for the creation of vortices in trapped 
BECs, although we considered only a few of them in this review. To illustrate the diversity of different methods, 
let us cite some other schemes. An experimental setup for vortex creation by Berry's phase induced Bose-Einstein 
condensation is proposed by Olshanii and Naraschcwski 17£]. A related vortex-production scheme employing the 
Aharonov-Casher effect is discussed by Petrosyan and You 180 1 . Other proposals suggest the creation of the vortex 
state by opto- mechanical stirring |167|; by a rotating force [ 181 1 ; by an adiabatic population transfer of a condensate 



from the ground to the excited Bose-condensed state via a Raman transition induced by laser light [166,168,165]; the 
accidental generation of vortices in a quench [ 182, 183 1 or in self-interference measurements [169]. 

A possible way to crea te ro tating states from a trapped ground-state BEC by using light-induced forces is proposed 
by Marzling and Zhang [184]. They show that the dipole potential induced by four traveling- wave laser beams with 
an appropriate configuration in space, phase and frequency can be used to realize such a system. Vortex states 
can be trappe d in an evaporative cooling process if the evaporation length is less than the size of the thermally 
excited state 185 185 1- In order to nucleate vortices, the trapped gas can be rotated at temperatures above the 



BEC transition. Rece ntly, it has been suggested that vorticity could be imprinted by imaging a BEC through an 
absorption plate [186]. The method consists of passing a far-off- resonant laser pulse through an absorption plate 
with an azimuthally dependent absorption coefficient, imaging the laser beam onto a BEC, and thus creating the 
corresponding nondissipative Stark-shift potential and condensate phase shift. A vortex ring may be formed by 
translating one condensate through another one [187] (this process is analogous to ring nucleation by moving ions 
in superfluid 4 He [188]), or by three-dimensional soliton decay [ 189, 190| . Recently the JILA group generated vortex 
rings by the decay of dark solitons through the snake instability 1 191]. 

Many different proposals for the detection of vortices in BECs have been mentioned in literature. Some of them 
are used in current experiments. The spatial size of the vortex core in the TF regime is too small to be observed; 
for visualizing the vortex state, it was suggested to switch off the trap and let the cloud expand ballistically fj2[ . 
After expansion, the size of the vortex core is magnified by approximately the same factor as the size of expanding 
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condensate ]7C,192|, so that the core becomes observable. Also the vortex state can be detected by the splitting of 
the collective condensate modes in axisymmetric traps [[32 111.115] or by looking at the phase slip in the interference 
fringes produced by two expanding condensates [167,70|. Dobrek et al. [ 186 1 proposed an interference method to 
detect vortices by coherently pushing part of the condensate with optically induced Bragg scattering. A detection 
scheme that reveals the existence of vortex states in a cylindrically symmetric trap is discussed by Goldstein, Wright 
and Meystre [ 193 1 . This scheme relies on the measurement of the second-order correlation function of the Schrodinger 
field and yields directly the topological charge of the vortex state. 

Also one can detect the vortex state by observing the off-resonance absorption image of the rotational cloud [ 166 1 . 
For a vortex state one should expect a bright "hole" in the image which accounts for the vortex core in the density 
distribution. Another possibility is to observe the Doppler frequency shift due to the quantized circular motion of the 
atoms [166], or by scattering fast atoms in a pure momentum state off a trapped atomic cloud [194|. 

Another question that has recently attracted significant theoretical and experimental interest is the dynamics and 
stability of dark soliton s an d vortex solitons in trapped condensates. Solitary waves (kin ks) h ave been studied in 
many physical contexts [195] and exist in different physical, c hemical a nd biological systems 196 . Recent theoretical 



studies discuss the dynamics and stability of dar k sol itons 197-201] (the range of parameters where the solitons 
are dynamically stable has been determined in f99],|l0l|| , while the theory of dissipative dynamics of a kink at finite 
temperature condensates has been developed in |202[| ), as well as suggestions for their cr eation 165,203,186]. Recently 
dark solitons inside a condensate were generated by a phase- imprinting method |204, 205[| . Unlike vortices, dark solitons 
are not topologically stable. At finite temperature, they exhibit thermodynamic and dynamic (small-amplitude) 
instabilities. The interaction of the soliton with the thermal cloud causes dissipation that accelerates the soliton. 
There is an interesting analogy between solitons and relativistic particles, in which the soliton velocity and speed of 
sound correspond to the particle velocity and speed of light in vacuum [202]. However, the kinematic mass of the 
soliton decreases when its velocity increases. This behavior is opposite to the case of relativistic particle, where the 
kinematic mass increases with velocity, and an infinite force is required to accelerate the particle beyond the velocity 
of light. In contrast to the particle, the soliton can reach the velocity of sound. An interesting problem is to create 
a soliton and a vortex simultaneously (this object is known as vortex soliton). The vortex soliton has a topological 
charge and therefore could be stable. 

Another challenging perspective for future experiments is the creation of vortex-like states in optically confined 
BECs. By relaxing the condition of spin polarization imposed by magnetic trapping, this new method of confinement 
permits the study of diverse textures that can be formed by the spinor order parameter, like those in superfluid 3 He-A 
p0| . Also optical traps allow strong variation of the scattering length via Feshbach resonances, which provides new 
possibilities for manipulating the condensate states. 

Among other challenging problems, one should mention measurements of vortex normal modes at higher temper- 
atures, which could establish the connection between the Bogoliubov approximation and self-consistent mean field 
theories. Also, it would be interesting to observe vortex dissipation and damping of vortex normal modes. 
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